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ATMOSPHERIC-ABSORPTION  ADJUSTMENT  PROCEDURE 
FOR  AIRCRAFT  FLYOVER  NOISE  MEASUREMENTS 


INTRODUCTION 

Working  Group  Sl-57  of  the  Standards  Committee  SI  on  Acoustics  developed 
a new  procedure  for  calculating  the  absorption  of  sound  by  the  atmosphere. 

In  1977,  the  new  procedure  was  in  the  process  of  being  approved  and  published 
by  the  American  National  Standards  Institute  (ANSI)  and  the  Acoustical  Society 
of  America  as  an  American  National  Standard  Method  for  the  Calculation  of  the 
Absorption  of  Sound  by  the  Atmosphere.  When  issued,  the  new  American  National 
Standard  will  have  the  ANSI  designation  SI. 26  and  hence  will  be  referred  to  in 
this  report  as  an  ANS  SI. 26.  The  technical  basis  and  the  experimental  verifi- 
cation for  the  new  calculation  procedure  are  documented  in  reports  prepared 
by  Shields  and  Bass  * and  by  Sutherland  ^ . 


The  new  procedure  for  calculating  atmospheric  absorption  is  based  on 
fundamental  physical  principles.  Extensive  laboratory  tests*  have  validated 
the  method  for  pure  tone  sounds  over  a wide  range  of  air  temperature  and 
humidity.  Although  the  primary  purpose  of  the  new  .ANSI  Standard  is  to  specify 
methods  for  determining  absorption  of  discrete-frequency  pure-tone  sounds  by 
the  atmosphere, the  proposed  Standard  also  presents  a general  description  of  the 
application  of  the  calculation  procedure  to  bands  of  noise,  e.g.,  sounds 
analyzed  by  constant-percentage-bandwidth  filters  such  as  1/ 3-octave-band 
filters . 

Techniques  applicable  to  each  particular  class  of  sound  propagation  problems 
are  not  necessarily  provided  in  ANS  SI. 26  because  the  new  Standard  is  intended 
to  be  basic  and  general  in  accordance  with  the  policy  of  the  Acoustical  Society 
of  America.  Thus,  it  may  be  necessary  for  a specific  user  to  develop  a method 
to  apply  the  proposed  procedures  to  a particular  problem,  such  as  analysis  of 
aircraft  flyover  noise  measurements. 


The  existing  procedure  for  adjusting  aircraft  flyover  noise  levels  for 
differences  in  atmospheric  absorption  losses  is  based  on  a combined 
theoretical,  experimental,  and  empirical  method  developed  by  the  Aircraft 
Noise  Measurement  Committee  A-21  of  the  Society  of  Automotive  Engineers  (SAE). 
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The  latest  version  of  the  SAE  procedure  is  contained  in  SAE  Aerospace 
Recommended  Practice  ARP  866A^ . The  procedures  of  SAE  ARP  866A  are  incorporat- 
ed directly  in  Part  36  of  the  Federal  Aviation  Regulations  (FAR  36)**,  in 
Annex  16  to  the  Standards  and  Recommended  Practices  of  the  International 
Civil  Aviation  Organization*,  and  in  International  Standard  IS  3891  of  the 
International  Standards  Organization  as  part  of  a procedure  for  describing 
noise  around  an  airport*. 

The  empirical  adjustments  to  laboratory  data  and  available  theory  that 
were  used  in  the  development  of  SAE  ARP  866  were  derived  from  aircraft  noise 
levels  and  meteorological  data  measured  a short  distance  above  ground  level. 

A basic  assumption  in  analysis  of  the  flyover  noise  test  results  was  that  the 
meteorological  measurements  near  the  ground  represented  conditions  all  along 
each  sound  path.  Based  on  a study  reported  by  McCollough  and  True^,  the  FAA ® 
has  proposed  modifying  the  original  procedures  of  FAR  36  to  require  applica- 
tion of  atmospheric  absorption  adjustments  by  the  SAE  ARP  866A  method  along 
segments  of  the  sound  propagation  paths.  The  FAA  'layered-weather'  or 
stratif ied-atmosphere  approach  is  intended  to  extend  the  use  of  SAE  ARP  866A 
to  varying  meteorological  conditions  aloft. 

This  report  describes  a computer  program  for  applying  the  proposed  ANSI 
Standard  ANS  SI. 26  to  the  problem  of  adjusting  measured  aircraft  flyover  noise 
levels  for  differences  in  atmospheric  absorption  losses  between  test  and 
reference  meteorological  conditions  along  the  propagation  path  connecting  the 
source  and  the  receiver.  The  computer  program  is  written  in  the  FORTRAN  IV 
language. 

The  following  Sections  (a)  describe  the  physical  characteristics  of  the 
problem  of  determining  adjustments  for  atmospheric  absorption  losses,  (b) 
give  reference  meteorological  conditions,  Ic)  describe  the  procedures  for 
determining  atmospheric  absorption  losses  by  the  methods  of  FAR  36  and  the 
proposed  ANS  SI. 26,  (d)  describe  the  development  of  the  computer  program,  (e) 
demonstrate  the  application  of  the  procedure  of  ANS  SI. 26  for  actual  aircraft 
noise  measurements,  (f)  outline  the  general  requirements  for  adapting  the 
computer  program  to  routine  processing  of  aircraft  noise  measurements,  and 
(h)  present  recommendations  for  investigations  that  should  be  conducted  prior 
to  considering  use  of  the  new  ANSI  procedures  for  aircraft  noise  type 
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certification. 


Appendices  contain  the  analytical  basis  for  the  adjustment  procedures 
for  discrete-frequency  pure-tone  sounds  and  for  bands  of  noise.  A listing 
of  the  computer  program  source  statements  and  definitions  of  all  symbols  are 
also  provided  in  Appendices. 

In  consonance  !,th  stated  Government  policy,  all  physical  quantities  are 
in  the  international  system  of  units.  Conversion  to  U.S.  customary  units,  if 
required,  may  be  made  by  reference  to  the  Metric  Conversion  Factors  given  in 
the  prefatory  material  to  this  report. 

DESCRIPTION  OF  THE  PROBLEM 

The  problem  of  determining  atmospheric  absorption  loss  adjustments  from 
test  to  reference  meteorological  conditions  is  illustrated  schematically  in 
figure  1.  Measured  sound  pressure  levels  are  shown  by  the  solid  line  in  the 
idealized  1/ 3-octave-band  spectrum  plot  on  the  right.  The  measured  levels 
represent  data  recorded  at  some  instant  of  time  as  the  airplane  flew  overhead 
with  airspeed  V and  Mach  number  M ; the  airplane's  position  on  the  flight 
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path  (shown  on  the  left  in  figure  1 as  a level-flight  flyover)  is  denoted  by 
the  angle  between  the  flight  path  and  a line  to  the  microphone.  The 
location  of  the  sound  source,  at  the  time  of  emission,  is  at  an  angle  ijy. 

Meteorological  conditions  at  the  time  of  the  flyover  noise  recording 
could  be  as  indicated  by  the  measured  profiles  of  temperature,  relative  , 

humidity,  and  pressure.  The  example  in  figure  1(a)  indicates  a strong  j 

ground-based  temperature  inversion  with  a weak  Inversion  aloft  at  a height  j 

above  the  height  of  the  aircraft.  The  humidity  profile  indicates  sub-  1 

stantially- lower  relative  humidity  aloft  than  at  ground  level  over  most  of  | 

the  sound  propagation  path.  The  low  upper-air  humidity  could  result  in 
significant  atmospheric  absorption  effects  on  sound  signals  propagating  from 
the  airplane  to  the  microphone.  The  air  pressure  measurements  show  the 
expected  trend  of  a reduction  in  pressure  with  Increasing  height.  Height  is 
assumed  to  be  measured  above  the  elevation  of  mean  sea  level.  i 

\ 

During  data  reduction,  the  recorded  sound  signals  are  analyzed  by  some 
type  of  filter  system  and  averaged  by  some  kind  of  detector  circuit.  Filter 
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systems  can  be  analog  or  digital.  As  indicated  by  the  sketch  in  figure  1(a), 

the  transmissibility  function  (filter  frequency  response)  of  practical  analog  | 

i 

filters  will  have  some  finite  transmission  in  the  stopbands  below  and  above 
the  cutoff  frequencies  at  the  limits  of  the  passband.  These  filter  skirts 
will  transmit  some  energy  and  increase  the  indicated  band  level  above  that 
which  would  be  indicated  by  an  ideal  filter. 

The  essence  of  the  problem  is  shown  in  figure  1(b)  — namely,  to  determine 
what  the  sound  pressure  levels  at  the  microphone  would  have  been  had  the  flyover 
noise  test  been  conducted  with  reference  rather  than  test  meteorological  con- 
ditions. Reference  meteorological  conditions  are  Indicated  by  defined  values 
for  temperature,  relative  humidity,  and  pressure  at  some  reference  height 
along  with  reference  lapse  rates.  The  reference  height  is  indicated  by  the 
square  symbols. 

The  frequency  response  of  the  filter  under  reference  conditions  is 
identical  to  the  frequency  response  under  the  test  meteorological  conditions. 

The  result  of  applying  atmospheric  absorption  corrections  generally 
yields  an  increase  in  high-frequency  sound  pressure  levels  as  indicated 
by  the  dashed  line  in  the  spectrum  sketch- 

The  balance  of  this  report  describes  development  of  a computer  program 
that  applies  the  procedure  of  the  proposed  ANS  SI. 26  to  adjust  aircraft 
flyover  noise  measurements.  The  specific  purpose  of  the  calculation  procedure 
is  to  adjust  1/3-octave-band  sound  pressure  levels,  measured  at  a point  near 
the  ground  during  an  aircraft  flyover  noise  test,  for  the  difference  between 
the  atmospheric  absorption  losses  that  occurred  at  the  time  of  the  test  and 
those  that  would  have  occurred  if  the  aircraft  noise  signal  had  propagated 
through  a reference  atmosphere  having  specified  temperature,  humidity,  and 
pressure  as  a function  of  height. 
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REFERENCE  METEOROLOGICAL  CONDITIONS 


FAR  36  Definitions 

As  In  all  earlier  versions,  the  October  1977  version  of  FAR  36**  defines 
reference  meteorological  conditions  in  paragraph  A36. 3(c)(1)  of  Appendix  A as 
follows 

(a)  sea  level  pressure  of  2116  psf  (76  cm  Hg)  (101.325  kPa) 

(b)  ambient  temperature  of  77°F  (25°C) 

(c)  relative  humidity  of  70  percent 

(d)  zero  wind. 

Appendix  A of  FAR  36  does  not  explicitly  state  a height  above  ground 
level  for  the  reference  pressure,  temperature,  and  relative  humidity.  Also, 
no  mention  is  made  of  corresponding  reference  lapse  rates  with  height  for 
pressure,  temperature,  relative  humidity,  or  wind. 

The  reference  meteorological  conditions  were  probably  meant  to  apply 
everywhere  along  the  sound  propagation  paths  because  the  original  development 
of  SAE  ARP  866  only  used  'surface*  weather  measurements  to  adjust  from  test 
to  reference  conditions.  Surface  weather  data  were  sometimes  measured  at  a 
height  above  ground  level  of  about  1.5  m somewhere  in  the  vicinity  of  the 
microphones,  and  were  sometimes  obtained  from  airport  facilities  where  the 
sensors  were  usually  about  10  m above  ground  level. 

No  specification  of  reference  lapse  rates  was  required  in  the  original 
version  of  FAR  36  because  atmospheric  absorption  adjustments  used  the  'surface' 
data  exclusively.  It  was  always  considered,  however,  that  the  'surface' 
meteorological  measurements,  wherever  measured,  should  be  representative  of 
the  conditions  near  the  microphones  at  the  time  of  the  aircraft  noise 
recordings.  Paragraph  A36.1(d)(4)  requires  that  the  airport  tower,  or  another 
facility,  must  be  approved  for  use  as  the  location  at  which  measurements 
are  made  of  atmospheric  parameters  representative  of  conditions  existing 
over  the  geographical  area  in  which  aircraft  noise  measurements  are  made. 
However,  the  surface  wind  velocity  and  temperature  must  be  measured  near  the 
microphone  at  the  approach,  sideline,  and  takeoff  measurement  locations. 


r " 

The  general  meteorological  test  conditions  defined  in  paragraph  A36.1(b)(3) 
for  acceptable  test  conditions  include  limits  on  ambient  temperature  and  wind 
at  a height  of  10  m above  ground  level.  No  height  is  specified  for  the 
acceptable  limits  on  relative  humidity.  It  has  become  general  practice,  with 
FAA  approval,  to  measure  'surface'  meteorological  test  conditions  of  pressure, 
temperature,  relative  humidity,  and  wind  at  a height  of  10  m. 

The  only  limitation  in  FAR  36  on  meteorological  conditions  aloft  is  the 
statement  in  A36.1(b) (3) (v)  that  there  should  be  no  temperature  inversion  or 
anomalous  wind  conditions  that  would  significantly  affect  the  noise  level  of 
the  aircraft  when  the  noise  is  recorded  at  the  measuring  points  defined  in 
Appendix  C.  The  wording  of  this  paragraph  was  the  basis  for  the  investiga- 
tions conducted  by  the  FAA  in  1974  and  reported  in  reference  7.  The  extensive 
changes  that  were  proposed  in  NPRM  76-21®  for  the  requirements  for  meteorolo- 
gical test  data  and  the  data-adjustment  procedures  also  stem  from  the  varying 
interpretations  that  were  permitted  by  A36 . 1(b) (3) (v) . 

Reference  Lapse  Rates 

In  order  to  calculate  differences  in  atmospheric  absorption  losses 
under  test  and  reference  conditions,  it  is  necessary  to  define  reference 
atmospheric  conditions.  Specification  of  reference  atmospheric  conditions 
along  the  sound  propagation  paths  requires  specification  of  reference 
lapse  rates  and  reference  conditions  at  some  reference  height. 

Lapse  rates  for  temperature  and  pressure  have  been  developed  and 
standardized.  Reference  9 defines  a standard  atmosphere  for  Northern 
midlatitudes  (45°  north)  and  average  annual  conditions.  References  10  and  11 
provide  atmospheric  data  at  other  latitudes  and  for  seasonal  variations. 

The  standard  atmosphere  of  reference  9 is  dry.  Relative  humidity  lapse 
rate  information  is  given  in  references  10  and  11. 

The  lapse  rate  for  temperature  is  -6.5  K/km  to  a height  above  mean  sea 
level  of  10  km.  The  relative  humidity  lapse  rate  is  approximately  -6.5 
percentage  points  per  km  to  a height  above  m.-an  sea  level  of  4 km.  For 
pressure  in  standard  atmospheres  (l.e.,  the  ..atic  of  the  air  pressure  P to 
the  standard  sea-level  atmospheric  pressure  Pq) » an  approximate  lapse  rate 
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for  the  logarithm  of  (P/P^),  derived  from  the  tabulated  values  In  reference  9, 

Is  -5.393*10  ^ log  (P/Pq)  per  meter  to  a height  above  mean  sea  level  of  4 km. 

To  be  consistent  with  the  definitions  of  FAR  36  for  reference  meteorological 
conditions  and  with  current  practice  for  a measurement  height  for  surface 
meteorological  data,  the  following  reference  conditions  were  defined  at  a 
reference  height  (altitude)  of  10  m above  mean  sea  level: 

• air  temperature:  298.15  kelvlns  (25.0°C); 

• relative  humidity:  70.0  percent;  and 

• air  pressure:  one  standard  atmosphere  or 
101.325  kPa  or  29.921  In.  Hg. 

These  definitions  for  reference  meteorological  conditions  at  the  10-m 
reference  height  are  not  the  same  as  those  in  references  9 and  11  where,  at 
the  reference  height  of  mean  sea  level,  the  standard  air  temperature  is 
288.15  kelvlns  (15°C),  the  standard  relative  humidity  is  75  percent,  and  the 
air  pressure  is  one  standard  atmosphere. 

Equations 


For  purposes  of  this  study,  reference  meteorological  conditions  were 
specified  by  the  following  equations: 


RTK  = RTKO  + RTKL  (H  - RH) ; 

(1) 

RHR  = RHRO  + RHRL  (H  - RH) ; 

(2) 

and 

RPA  = 10^^^ 

(3) 

where 

LPA  = RPAL  (H  - RH). 

(4) 

The  symbols  have  the  following  meanings.  RTK  Is  the  reference  temperature. 
In  kelvlns,  and  RHR  Is  the  reference  relative  humidity,  in  percent,  at  any 
height  H.  RTKO  and  RHRO  are  the  reference  values  at  the  reference  height 
(i.e.,  298.15  K and  70  percent).  RTKL  and  RHRL  are  the  reference  lapse  rates 
(i.e.,  -0.0065  K/m  and  -0.0065  percentage  points  per  meter).  H is  height.  In 
meters,  above  mean  sea  level;  RH  Is  the  reference  height  of  10  m.  RPA  Is  the 
reference  air  pressure  in  standard  atmospheres  at  any  height.  LPA  is  the  log 
of  the  pressure  in  standard  atmospheres  at  any  height  and  RPAL  is  its  lapse 
rate  of  -5.393  x 10  ^ per  meter. 
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Note  that  uniform  meteorological  conditions  could,  if  desired,  be  generat- 
ed by  setting  RTKL,  RHRL,  and  RPAL  equal  to  zero. 

s FAR  36  ATMOSPHERIC-ABSORPTION  ADJUSTMENT  PROCEDURE 

1 

! Paragraph  A36. 3(d)(1)  of  FAR  36  requires  that  measured  noise  data  be 

adjusted  from  the  test  meteorological  conditions  [measured  in  accordance  with 
r . the  requirements  of  paragraph  A36.1 (d)(4)]  to  the  reference  meteorological 

f conditions  of  paragraph  A36.3(c).  The  adjustments  must  be  made  using  the 

methods  of  the  latest  version  of  SAE  ARP  866  by  the  general  procedures  given 
in  paragraph  A36.5.  Detailed  data-correctlon  procedures  are  given  in  para- 
I , graph  A36.6. 

I The  detailed  procedures  of  A36.6  determine  a correction  factor  to  be 

added  to  the  measured  value  of  the  test-day  effective  perceived  noise  level 
(EPNL) . The  correction  factor  is  intended  to  account  (1)  for  the  change  in 
atmospheric  sound  absorption  over  the  distance  from  the  airplane  flight  path 
to  the  microphone  at  the  time  associated  with  the  maximum  value  of  the  tone- 
corrected  perceived  noise  level  (i.e.,  at  the  time  of  occurrence  of  the 
test-day  PNLTM) , and  (2)  for  the  change  in  atmospheric  sound  absorption  caused 
by  differences  in  the  sound  propagation  pathlength  between  the  location  of  the 
microphone  and  the  locations  of  the  effective  sound  source  on  the  test  and  the 
reference  airplane  flight  paths. 

The  correction  factor  that  is  defined  in  paragraph  A36.6(a)  to  account 
for  differences  between  test-day  and  reference-day  atmospheric  absorption  is 
A1  = PNLT^gf  - PNLTMj.ggj.  = AEPNL  = EPNLj.gf  - EPNL^gg^.  The  value  of  tone- 
corrected  perceived  noise  level  under  acoustic  reference-day  conditions 
(PNLTj-gf)  is  determined  by  adjusting  the  measured  1/ 3-octave-band  sound 
pressure  level  spectrum  associated  with  PNLTMteg^.  The  adjustment  procedure 
> accounts  for  differences  in  atmospheric  absorption  and  acoustic  pathlength  by 

‘ means  of  separate  formulas  given  in  paragraph  A36.6(d)  for  the  takeoff,  sideline,  j 

I and  approach  noise  measurements.  j 

I ! 

I Changes  in  the  duration  correction  factor  are  determined  from  the  ratio  j 

of  the  minimum  distances  from  the  microphone  to  the  test  and  the  reference 
flight  paths. 

9 


Microphones  are  required  by  FAR  36  to  be  placed  at  a height  of  1.2  m 
above  the  ground  surface.  Interference  effects  between  sound  waves  propagating 
directly  from  the  source  to  the  microphone  and  sound  waves  reflected  from  the 
ground  surface  introduce  peaks  and  nulls  in  the  spectra  at  frequencies  that 
depend  on  the  geometry  of  the  test  setup  and  the  speed  of  the  source.  The 
magnitude  of  the  spectral  irregularities  caused  by  ground  reflection  effects 
can  be  several  decibels.  The  usual  industry  practice  is  to  not  remove 
irregularities  caused  by  ground  reflections  before  calculating  atmospheric 
absorption  adjustments  to  the  measured  noise  levels. 

The  spectrum  of  the  aircraft  noise  signal  at  the  microphone,  compared 
with  the  spectrum  of  the  noise  at  the  source,  is  shifted  to  higher  or  lower 
frequencies  by  an  amount  that  depends  on  the  airplane  Mach  number  M and  the 
sound  propagation  angle  i//.  This  Doppler  effect  causes  an  apparent  shift  of 
the  source  spectrum  to  higher  frequencies  as  the  airplane  approaches  the 
microphone  [by  1/(1-M  cos  yj;)  ] and  to  lower  frequencies  after  the  airplane 
passes  the  microphone  [by  1/(1+M  cos  ij^)  ] . For  the  range  of  airplane  Mach 
numbers  and  propagation  angles  normally  encountered  in  flyover  noise  testing, 
the  frequency  shift  is  on  the  order  of  10  to  25  percent. 

For  broadband  random  sounds,  it  is  not  necessary  to  consider  Doppler 
frequency  shifts  to  the  observed  spectra.  For  discrete-frequency  signals, 
however,  one  should  consider  the  effect  of  the  source  motion  when  attempting 
to  relate  the  frequency  of  the  source  (e.g.,  a fundamental  fan-blade-passage 
frequency)  to  the  frequency  of  a signal  identified  as  a pure  tone  in  the 
passband  of  some  1/ 3-octave-band  filter.  Knowledge  of  the  true  frequency  of 
the  pure-tone  source  is  needed  to  determine  the  pure-tone  atmospheric- 
absorption  adjustment  factor  at  the  Doppler  shifted  frequency. 

ATMOSPHERIC-ABSORPTION  ADJUSTMENT  PROCEDURE  BY  ANS  SI. 26 

Each  0.5-s  sample  of  the  aircraft  noise  signal  is  assumed  to  be  separated, 
if  required,  into  broadband  and  discrete-frequency  spectral  components.  An 
algorithm  is  assumed  to  be  available  to  search  the  measured  spectra  and  to 
determine  the  separate  components.  Atmospheric-absorption  adjustments  are 
determined  separately  for  each  spectral  component.  Details  of  the  analytical 
basis  for  the  pure-tone  and  the  broadband  adjustment  procedures  are  given  in 
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Appendices  A and  B,  respectively.  The  discussion  here  gives  only  the  general 
outline  of  the  methods. 

Pure  Tones 

The  propagation  path  is  assumed  to  be  divided  into  a number  of  segments 
over  which  an  average  temperature,  humidity,  and  pressure  can  be  established. 
The  atmospheric  absorption  at  the  frequency  of  the  tone  (i.e.,  at  the  Doppler- 
shifted  frequency)  is  calculated  first  for  the  meteorological  conditions 
of  the  atmosphere  at  the  time  of  the  test  and  then  for  the  meteorological 
conditions  of  the  reference  atmosphere.  The  difference  is  added  to  the 
measured  pure-tone  sound  pressure  levels  to  obtain  the  sound  pressure  levels 
adjusted  to  reference  meteorological  conditions. 

The  pure-tone  adjustment  factor,  AL  „„  , is  defined  as 

1 9 XK 


^4,TR  = ^ 


where  the  adjustment  factor  A is  determined  from 


A = 10 


[(a-  .-a  )(AC,/10)  + • • • + (a  -a  )(AC  /lO)] 

1,11,2  1 n,l  n,2  n 


The  lengths  of  the  various  segments  of  the  propagation  path  are 
AC2,  *••  to  The  pure-tone  absorption  coefficients  a for  the  different 

meteorological  conditions  on  each  path  segment  are  calculated  by  the  method 
of  the  proposed  ANS  SI. 26  at  the  specified  frequency.  Equations  (A2)  to  (A6) 
from  Appendix  A are  used  to  determine  the  absorption  coefficient  a in  nepers 
per  meter.  The  absorption  coefficient  a in  decibels  per  meter  is  then 
determined  from  a = 8.686  a.  The  lengths  of  the  propagation  path  segments 
are  determined  by  dividing  the  length  of  the  total  sound  propagation  path  into 
Intervals  according  to  the  heights  at  which  meteorological  data  are  measured. 

Bands  of  Noise 

Development  of  the  analytical  basis  for  determining  atmospheric  absorp- 
tion losses  for  bands  of  wideband  random  noise  and  the  corresponding  computa- 
tional procedure  were  the  main  technical  efforts  of  the  study  reported  here. 
The  general  procedure  involves  the  evaluation  of  two  integrals  over  the 
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frequency  range  of  the  passband  of  each  of  the  contiguous  filter  bands. 


The  adjustment  factor,  AL  for  a band  of  noise  at  center  frequency 

D,  TK  1 

fj^  with  passband  lower  and  upper  frequencies  f^^  and  fy  is  determined  from 

ALg^TR^^i^  = 10  log  {[^f  a df]/[^/  df]}  (7) 

li  L 

where  A is  the  absorption  factor  of  equation  (6)  over  the  propagation  path  and 
Gj^  is  the  spectral  density  of  the  sound  pressure  signal  at  the  microphone 
before  passage  through  the  data-acquisitlon  and  data-processlng  system. 

Equation  (7)  assumes  that  the  frequency  response  of  the  real  filters  can  be 
adequately  represented  by  the  response  of  ideal  filters. 

The  central  problem  in  evaluating  equation  (7)  is  how  to  determine  a 
reasonable  approximation  to  G^^  knowing  only  the  measured  sound  pressure  levels 
after  passage  through  the  data-acquisitlon  and  data-processing  system.  For 
each  passband,  the  approximation  chosen  for  the  analysis  in  Appendix  B is 

G^(f)  = G^(f^)  (f/f.)^'  (8) 

where  is  the  pressure  spectral  density  at  the  band  center  frequency  and 

£'  is  the  slope  of  the  spectrum.  On  a plot  of  log  G^^  vs  log  f,  the  approxima- 
tion of  equation  (8)  is  a straight  line  with  slope  I'  over  the  frequency  range 
from  fy  to  fy. 

With  equation  (8),  equation  (7)  becomes 

ALg,TR(fi)  = 10  log  A f^’  df]/[^/  ^ f^'  df]}.  (9) 

L li 

The  numerator  term  in  equation  (9)  is  evaluated  as  the  summation  of  a 
series  of  consecutive  integrals.  The  number  of  Integrals  in  the  sum  is 
determined  by  dividing  the  filter  passband  into  SE  segments  or  constant- 
percentage-bandwidth  frequency  Intervals. 

To  determine  an  Integrable  expression  for  A,  it  is  assumed  that  the 
absorption,  over  any  of  the  segments  AC  which  make  up  the  sound  propagation 
path,  is  a function  of  frequency  only.  Over  the  narrow  range  of  frequency 
of  one  of  the  Integrals  in  the  sum,  the  absorption  is  assumed  to  be 


12 


proportional  to  some  power  of  the  frequency  f.  This  assumption  leads  to  the 
expression 


A(f)  = A(fj)  [(f/f^)*"] 


where  ^(fj)  is  the  pathlength  absorption  factor  at  the  frequency  f^  at  the 
beginning  of  the  frequency  range  for  one  of  the  integrals  and  K is  an  effective 
absorption- spectrum  slope  parameter  averaged  over  the  pathlength  and  over  the 
differences  between  absorption-spectrum  slopes  under  test  and  reference 
meteorological  conditions. 

The  slope  i'  of  the  sound  pressure  spectral  density  is  estimated  from  the 
band-level  noise-spectrum  slope  H using 

i'  = i - 1.  (11) 

With  equations  (10)  and  (11),  evaluation  of  equation  (9)  leads  to  four 
possible  cases. 

When  1^0,  the  band-loss  adjustment  factor  is 


where 


"S,TR^^i>  “ ih 


BLJl  = BLlj  when  K+i?i0 


BLJl  = BL2j  when  K+)l=0. 


When  1=0  and  K?^0,  then 


When  K=0,  then  there  is  no  difference  between  the  test  and  reference 


conditions  and 


regardless  of  noise  slope  H. 


Appendix  B gives  the  specific  expressions  for  the  terms  j » and 

BL3j  as  well  as  specific  procedures  for  determining  the  band  center  frequencies. 
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the  passband  limiting  frequencies,  and  a choice  for  the  number  of  passband 
frequency  intervals  SE.  Specific  recommendations  are  also  given  for 
determining  the  noise  slope  i for  each  band. 

General  Limitations 

In  agreement  with  the  recommendations  of  the  proposed  ANS  SI. 26,  applic-  ^ 

ability  of  the  procedures  described  here  for  pure  tones  or  bands  of  noise  is 
limited  to  air  temperature  in  the  range  from  0.0°  to  40.0°  Celsius,  relative 
humidities  from  10.0  to  100.0  percent,  and  nominal  1/ 3-octave-band  center 
frequencies  from  50  Hz  to  10  kHz.  The  procedures  are  also  limited  to  air  of 
standard  composition  and  normally  encountered  air  pressures. 

APPLICATION  TO  ADJUSTMENT  OF  AIRCRAFT  FLYOVER  NOISE  MEASUREMENTS  i 

1 

The  absorption-loss  adjustment  factors  can  be  applied  to  measured  spectra  ■ 

associated  with  any  data  sample.  The  procedure  is  thus  applicable  to  the  j 

FAR  36  method  of  adjusting  only  the  spectrum  associated  with  the  maximum  value 

i 

of  the  test-day  tone-corrected  perceived  noise  level  as  well  as  to  adjustment  ■ 

1 

of  every  0.5-s  data  sample  throughout  the  duration  of  a flyover  — as  might  be 
desired  for  research  investigations  of  aircraft  noise  propagation. 

1 

I 

General  Assumptions  | 

The  following  general  assumptions  were  made  in  the  calculation  of  absorp-  | 

tion-loss  adjustment  factors.  j 

S 

1.  For  each  frequency  of  interest,  the  effective  source  of  sound  from 
the  airplane  is  at  a single  point  on  the  flight  path  that  can  be  determined, 
for  each  instant  of  time  during  a flyover,  from  the  airplane’s  position  and 
airspeed.  The  airplane's  position  is  defined  by  the  coordinates  of  an 
airplane  reference  point.  The  coordinates  include  position  along  the  flight 
track,  lateral  deviation,  and  height  above  ground  level  (height  AGL) . 

2.  For  each  data  sample  throughout  the  flyover  noise  recording,  the 
direction  of  sound  propagation  from  the  source  to  the  receiver  (microphone) 
is  constant  and  defined  by  a propagation  angle  whose  sides  are  the  flight 
path  and  a straight  line  (sound  ray)  between  the  source  and  the  receiver. 

Bending  of  the  sound  rays  caused  by  refraction  effects  occurring  during 
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propagation  through  regions  or  layers  of  different  temperature  (different 
speeds  of  sound)  Is  negligible.  Bending  of  the  sound  rays  caused  by  propaga- 
tion through  winds  of  varying  speeds  and  direction  Is  also  negligible  because 
the  wind  Is  calm  at  all  points  along  every  sound  propagation  path. 

3.  Attenuation  and  scattering  effects  caused  by  sound  propagation 
through  atmospheric  turbulence  are  negligible  for  each  sound  ray. 

4.  Aircraft  noise  signals  are  filtered  Into  twenty  four  1/3-octave  bands 
with  nominal  band  center  frequencies  ranging  from  50  to  10,000  Hz. 

5.  Nonlinear  propagation  effects  caused  by  high  sound  pressure  levels 
at  the  source  are  negligible. 

6.  No  atmospheric  absorption  adjustments  are  to  be  made  when  the  aircraft 
noise  signal  is  not  5 dB  or  more  above  the  level  of  the  corresponding  ambient 
noise  in  any  1/3-octave  band. 


Assumptions  Regarding  Input  Data 

The  following  assumptions  were  made  for  the  data  required  as  input  to  the 
calculations.  The  assumptions  apply  to  the  general  problem  of  adjusting 
measured  aircraft  noise  signals  for  FAR  36  or  research  purposes  as  well  as  the 
particular  test  cases  selected  for  demonstration. 

Meteorological  Parameters.  — Temperature  and  humidity  data  are  measured 
at  regular  Intervals  of  height  AGL  at  several  times  throughout  each  test  day. 
Values  of  temperature  and  humidity,  applicable  to  each  set  of  test  sound 
pressure  levels  from  a flyover  noise  measurement,  can  be  interpolated  from  the 
available  measurements.  The  unit  for  temperature  is  kelvins;  for  relative 
humidity  it  is  percent.  Temperature  and  relative  humidity  profiles  are  avail- 
able to  a height  greater  than  the  airplane  height.  Height  is  measured  In 
meters. 


Station  barometric  pressure  at  a height  c^  10  m is  measured  periodically 
throughout  each  test  day.  The  pressure  Is  applicable  to  the  actual  runway 
elevation  and  not  corrected  to  sea  level.  A value  for  barometric  pressure 
applicable  to  the  time  of  the  flyover  noise  test  can  be  interpolated  from  the 
available  measurements.  The  unit  for  barometric  pressure  measurements  is 
the  kilopascal  (or  inches  of  mercury). 
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(If  air  pressure  measurements  as  a function  of  height  AGL  are  available, 
it  is  preferable  to  use  the  measured  pressure-profile  data,  interpolated  to 
the  time  of  the  flyover  noise  test,  in  conjunction  with  the  measured  tempera- 
ture and  relative  humidity  profiles  rather  than  the  single  10-m  pressure 
measurement . ) 

Airplane  Tracking.  — The  coordinates  of  an  airplane  reference  point  are 
available  as  a function  of  time  throughout  the  duration  of,  and  in  synchroni- 
zation with,  the  recording  of  the  airplane's  noise  signal.  Airplane  coordinates 
are  relative  to  an  origin  point  on  the  ground  and  consist  of  x,  y,  z data. 

The  x-coordinate  is  in  the  ground  plane  and  along  the  nominal  flight  track. 

The  z-coordinate  is  height  AGL.  The  y-coordinate  is  in  the  ground  plane  and 
is  a measure  of  lateral  deviation  from  the  nominal  flight  track.  The  unit  for 
distance  measurements  is  the  meter. 

Noise  Source.  — The  parameters  that  determine  the  acoustical  characteristics 
of  the  airplane  noise  source  are  constant  throughout  the  duration  of  each 
sampling  interval,  i.e.,  engine  power  setting,  airspeed,  airplane  configuration, 
and  attitude  are  all  constant  over  each  time  sample. 

Sound  Pressure  Levels.  — Sound  pressure  levels  in  1/3-octave  bands  are 
available  at  0.5-s  intervals  throughout  the  duration  of  each  flyover  noise 
recording.  Sound  pressure  levels  represent  true  mean-square  sound  pressures, 
are  in  decibels  (dB)  referenced  to  20  micropascals,  and  are  corrected  for  all 
known  sources  of  error  including:  windscreen  corrections,  microphone  response 
corrections,  recording  and  processing-system  frequency-response  corrections, 
system  sensitivity,  and  interference  by  ambient  noise. 

Ambient  Noise  Corrections.  — Corrections  for  ambient  noise  interference 
are  made  when  the  aircraft  noise  signal  is  5 dB  or  more  above  the  corresponding 
ambient  noise  signal  in  each  1/ 3-octave  band.  When  the  aircraft  noise  signal 
is  not  more  than  5 dB  above  the  ambient  noise  signal,  the  data  processing 
system  does  not  attempt  to  make  a correction  and  sets  the  1/3-octave-band 
sound  pressure  level  to  0.0  dB  or  some  other  negligible  value. 

Averaging  Time.  — other  than  integration  over  the  0.5-s  duration  of  the 
sampling  interval,  no  additional  smoothing  of  the  sound  pressure  level 
measurements  is  included  in  the  data,  e.g.,  no  running  averaging  of  the  data 
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samples  is  made  to  simulate  sound  pressure  levels  that  would  have  been  measur- 
ed on  a sound  level  meter  with  high  damping. 

Identification  of  Data  Sample  Time.  — The  instant  of  time  associated  with 
the  flyover  noise  data  samples  is  the  midpoint  of  the  0.5-s  sample  period 
because  that  time  is  most  representative  of  the  aircraft's  position  during 
the  0.5-s  period. 

The  assumptions  described  above  for  averaging  time  and  identification  of 
the  data  sample  time  may  not  comply  with  the  data-processing  techniques  used 
by  all  applicants  for  aircraft  noise  type  certification  under  FAR  36.  Similarly, 
not  every  applicant  may  use  the  5-dB  rule  for  rejecting  data  samples  contaminat- 
ed by  ambient  noise;  some  may  use,  for  example,  a 3-dB  rule.  The  specific  and 
the  general  assumptions  made  here,  however,  are  considered  to  be  a consistent 
interpretation  of  the  proper  application  of  atmospheric  absorption  corrections. 

Special  Considerations 

There  are  three  Important  considerations  that  should  be  taken  into  account 
in  determining  atmospheric  absorption  loss  adjustments  to  aircraft  flyover 
noise  measurements.  Development  of  techniques  to  overcome  the  problems  arising 
from  the  considerations  in  actual  applications  was  beyond  the  scope  of  the 
project . 

Ground  Reflection  Effects.  — Absorption  loss  adjustments  should  be  applied 
to  sound  spectra  that  represent  equivalent  free-field  sound  pressure  levels. 
Measurements  from  microphones  located  above  a ground  plane  will  have  ground- 
reflection  effects  caused  by  interference  between  direct  and  reflected  waves. 
Unless  the  microphone  height  is  large  with  respect  to  the  longest  wavelength 
of  interest  or  unless  the  microphone  is  mounted  flush  in  an  acoustically  large 
and  rigid  surface,  the  measured  spectra  will  have  Interference  peaks  and  nulls 
in  the  frequency  range  of  interest.  Spectral  irregularities  should  be  removed 
before  attempting  to  apply  atmospheric  absorption  corrections  because  the  methods 
of  the  proposed  ANS  SI. 26  cannot  be  used  to  account  for  ground  attenuation 
effects  on  reflected  waves. 

If  the  acoustical  impedance  of  the  ground  surface  is  known  (or  can  be 
reasonably  estimated),  then  available  analytical  techniques  can  be  applied  to 
remove  the  ground  reflection  effects.  If  significant  spectral  Irregularities 


17 


from  ground-reflection  effects  are  present  in  the  measured  data,  it  will  be 
necessary  to  remove  them  in  order  to  avoid  wrongly  Identifying  an  interference 
peak  as  a discrete-frequency  sound  from  the  airplane. 


Separation  of  Discrete  and  Broadband  Components.  — After  ground-reflection 
effects  have  been  removed  from  the  measured  sound  pressure  levels  to  yield  a 
spectrum  representative  of  free-field  sound  pressure  levels,  the  spectra  should 
then  be  examined  to  identify  the  frequencies  and  levels  of  discrete-frequency 
components.  Pure-tone  absorption  loss  corrections  are  applied  to  the  discrete- 
frequency  components;  band- loss  corrections  are  applied  to  the  filtered  broad- 
band levels.  The  tone  identification  procedure  of  reference  4 could  be  used  to 
establish  the  levels  of  the  separata  discrete  and  broadband  spectral  components. 

Missing  Band  Levels.  — To  avoid  attempting  to  apply  atmospheric  absorption 
corrections  to  invalid  data,  it  is  necessary  for  the  data  processing  system  to 
identify  the  bands  irretrievably  contaminated  by  ambient  noise  and  to  eliminate 
them  from  determination  of  noise  slopes  or  absorption  losses.  If  band  levels 
are  missing  from  only  the  beginning  or  end  of  the  set  of  24  bands  with  center 
frequencies  between  50  and  10,000  Hz,  then  the  calculation  of  absorption  losses 
could  proceed  by  limiting  the  analysis  to  those  bands  for  which  valid  data 
exist.  If  data  are  not  available  for  bands  scattered  between  50  and  10,000  Hz 
or  if  there  are  two  or  more  adjacent  (but  not  initial  or  final)  bands  with 
missing  levels,  then  either  special  rules  based  on  additional  assumptions 
would  need  to  be  invoked  (because  of  the  requirement  to  have  a value  for  the 
slope  of  the  noise  spectrum  as  a parameter  in  computing  the  band  losses)  or 
no  attempt  should  be  made  to  determine  absorption  losses  in  these  cases. 

DESCRIPTION  OF  COMPUTER  PROGRAM 

Adjustment  of  measured  aircraft  noise  levels  for  differences  between  test 
and  reference  meteorological  conditions  is  based  on  the  analytical  approach 
described  above  and  in  Appendices  A and  B.  For  any  sample  of  data  throughout 
the  duration  of  a measured  aircraft  noise  signal,  the  calculation  proceeds  by 
determining  the  atmospheric  absorption  loss  along  the  calculated  propagation 
path  from  the  microphone  to  the  source  of  sound  using  atmospheric  conditions 
existing  at  the  time  of  the  test  and  then  from  the  source  to  the  microphone 
using  reference  meteorological  conditions.  The  difference  between  the 
atmospheric  absorption  losses  along  the  path  determines  the  magnitude  and  sign 
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of  tiie  adjustment  required  for  each  1/ 3-octave-band  sound  pressure  level.  The 
computer  program  that  was  prepared  to  calculate  adjustment  factors  for  bands 
of  noise  is  called  BANDLOSS  and  is  described  in  Appendices  C and  D. 


Figure  2 indicates  the  steps  involved  in  implementing  a method  to  calcu- 
late band-loss  atmospheric  absorption  adjustment  factors  by  means  of  three 
sequentially  stepped  DO  loops  and  one  subroutine.  The  DO  loops  consist  of 
one  main  loop  and  two  inner  or  nested  loops. 

After  reading  in,  calculating,  or  defining  the  necessary  input  informa- 
tion, the  program  enters  the  outer  calculation  loop.  The  outer  loop  steps 
through  the  set  of  24  bands  comprising  the  measured  input  sound  pressure  levels. 
The  first  inner  loop,  loop  2,  performs  the  numerical  Integration  on  frequency 
fj  over  the  frequency  range  of  the  passband  of  each  filter  and  returns  the 
value  of  BLJl  or  BL3J  to  the  outer  loop. 

At  each  increment  of  f^  from  j = 1 to  j = SE,  the  first  inner  loop  calls 
the  second  inner  loop,  loop  3,  for  a calculation  of  the  effective  absorption- 
slope  parameter  K and  the  exponent  AE  of  the  absorption  factor  A(f=fj)  = AJ. 

The  third  loop  calculates  the  quantities  AE  and  K for  the  test  and  reference 


INPUT  INFORMATION 


Figure  2.-lllustration  of  the  calculation  of  atmospheric  absorption  adjustment 
factors  and  adjusted  1 /3-octave-band  sound  pressure  levels. 
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meteorological  conditions  applicable  to  each  segment  of  the  sound  propagation 
paths,  i.e.,  at  increments  in  height  from  H = 1 to  H = lA  where  lA  is  an  index 
counter  representing  the  height  for  the  meteorological  measurements  that  is 
just  less  than  the  height  of  the  aircraft  at  the  time  when  the  sound  was 
emitted  from  the  aircraft  noise  source.  At  each  step  in  height,  loop  3 calls 
subroutine  ANSAB  to  calculate  the  sound  absorption  coefficient  for  specified 
frequency,  temperature,  humidity,  and  pressure. 

The  total  number  of  calculations  required  by  the  three  loops  is  the  pro- 
duct of  the  ranges  of  the  DO  loops  or  24  x SE  x lA.  For  typical  aircraft 
noise  measurements,  the  index  lA  could  have  a value  ranging  from  15  to  30. 

The  number  of  frequency  segments  is  SE  = 20.  Thus,  the  three  loops  require  up 
to  15,000  calculations.  Subroutine  ANSAB  is  called  four  times  at  each  height 
step  by  loop  3. 

The  BUI  or  BL3J  sum  for  each  band  is  passed  back  to  loop  1 by  loop  2. 

The  band  loss  adjustment  factor  AL„  is  the  output  of  loop  1.  Adjusted 

D y IK 

sound  pressure  levels  are  obtained  by  adding  the  calculated  ALg  values  to 
the  input  sound  pressure  levels. 

Input  information  consists  of  the  various  quantities  outlined  in  figure  3. 
The  measured  inputs  are  the  noise,  meteorological,  and  airplane  data.  Calculat- 
ed input  data  consist  of  the  heights  for  the  meteorological  measurements,  the 
height  at  the  midpoint  of  each  height  interval  and  the  corresponding  air 
temperature,  relative  humidity,  and  atmospheric  pressure  in  standard  atmospheres. 
Temperature,  relative  humidity,  and  pressure  under  reference  conditions  are 
also  calculated.  Sound  propagation  path  lengths  and  the  differences  in  the 
levels  of  the  input  1/3-octave-band  sound  pressure  levels  are  determined. 

The  precise  values  of  the  band  center  frequencies  are  calculated. 

The  reference  meteorological  conditions,  the  band  frequency  ratio,  the 
number  of  frequency- interval  segments  in  a filter  passband,  and  the  International 
Standard  Band  Numbers  for  use  in  calculating  the  band  center  frequencies  are 
all  defined. 

Appendix  C lists  the  FORTRAN  source  statements  for  program  BANDLOSS,  in- 
cluding subroutine  ANSAB,  for  computation  of  atmospheric  absorption  adjustment 
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DEFINITIONS 

REFERENCE  METEOROLOGICAL  PARAMETERS:  RHM,  RTKO,  RTKL, 

RHRO,  RHRL,  RPAL, 
RBPO 

FREQUENCY  RATIO:  RF 

NUMBER  OF  FREQUENCY-INTERVAL  SEGMENTS:  SE 
INTERNATIONAL  STANDARD  BAND  NUMBERS:  ISBN 


Figure  3.-lnput  information  needed  for  calculation  of  ALg  jp  and  SPLA. 

factors  in  accordance  with  the  analyses  in  Appendices  A and  B.  Appendix  D 
describes  and  defines  the  symbols  and  terms  used  in  program  BANDLOSS. 

The  statement  line  numbers  in  Appendix  C for  subroutine  ANSAB  are 
sequential  to  those  of  the  main  part  of  program  BANDLOSS  because  the  subroutine 
and  the  main  program  were  merged  together  for  calculating  the  test  cases.  The 
statements  for  the  input  data  were  developed  specifically  for  the  demonstration 
test  cases  described  in  the  next  section.  Relatively  minor  modifications  would 
be  required  to  adapt  the  procedures  described  here  to  a specific  flyover-noise 
data-acqulsition/data-reductlon  system.  The  main  calculation  loops  could  be 
used  without  alteration.  The  program  was  written  for  direct  on-line  submittal 
from  an  Interactive  terminal  to  a large-scale  digital  computer.  No  major 
changes  would  be  required  for  batch-mode  submittal  through  a terminal,  or  for 
conversion  to  input  through  a card  reader,  or  for  incorporation  in  a data- 
processing  program. 
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DEMONSTRATION  OF  APPLICATION  OF  COMPUTER  PROGRAM 


To  demonstrate  application  of  the  computation  procedures  to  aircraft  noise 
test  data,  we  consider  the  situation  where  the  following  assumptions  are  valid. 

1.  The  airplane  flew  a level  flight  path  in  a straight  line  over  the 
microphone  (l.e.,  y = 0.0)  such  that  the  airplane  height  was  constant  throughout 
the  duration  of  each  noise  recording  (i.e.,  z = constant).  The  airplane  posi- 
tion was  synchronized  with  the  noise  recording  and  the  time  when  the  airplane 
was  directly  over  (or  closest  to)  the  microphone  is  known. 

2.  The  airplane  flew  over  a runway  whose  elevation  was  approximately  equal 
to  that  of  mean  sea  level. 

3.  The  spectrum  of  the  aircraft  noise  signal  for  each  time  sample  is 
equivalent  to  a free-field  spectrum.  No  corrections  are  required  for  spectral 
irregularities  caused  by  ground  reflection  effects.  The  calculation  procedure 
does  not  need  to  contain  a method  to  search  the  spectrum  and  identify  and  re- 
move spectral  irregularities  caused  by  ground  reflections. 

A.  The  noise  spectra  used  for  demonstration  of  the  atmospheric-absorp- 
tion correction  method  represent  broadband  sound  sources  only.  The  calcula- 
tion procedure  does  not  need  to  contain  a method  to  search  the  spectrum  and 
identify  the  frequency  and  level  of  pure-tone  components  and  the  remaining 
broadband  components . 

5.  The  demonstration  spectra  are  not  contaminated  by  interference  from 
ambient  noise.  Valid  sound  pressure  levels  are  available  for  all  2A  bands 
with  nominal  center  frequencies  from  50  to  10,000  Hz.  The  calculation 
procedure  does  not  need  to  contain  a method  to  search  the  spectrum  to  identify, 
count,  and  locate  bands  where  the  corresponding  ambient  noise  levels  were  too 
high  and  the  aircraft  noise  signal  was  set  to  the  negligible  value. 

6.  There  were  31  heights  (IMAX  = 31)  at  which  air  temperature  and 
relative  humidities  were  measured.  Air  pressure  was  measured  only  at  a height 
of  10  m. 

With  these  assumptions  the  geometrical  relationship  between  the  airplane 
and  the  microphone  is  as  shown  in  figure  A.  The  airplane  is  at  a height  AH 
above  ground  level;  the  microphone  is  at  a height  ZM.  The  height  between  the 
flight  path  and  the  microphone  is  AMH  “AH  - ZM. 
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Figure  4.-Geoinetrical  relationships  for  test  cases. 

When  the  airplane  is  at  a location  on  the  flight  path  defined  by  angle 

()>  to  the  microphone,  the  sound  was  emitted  from  an  effective  source  with 

emission  angle  ii  between  the  flight  path  and  the  ray  between  the  source  and 

the  microphone.  The  time  t_  when  the  sound  was  received  at  the  microphone 

is  measured  relative  to  the  time  when  the  aircraft  was  over,  or  closest  to, 

the  microphone,  i.e.,  where  4>  = 0 and  ij;  » cos  M with  M the  airplane  Mach 

d 3l 

number  (M  = V /c  where  V is  the  airspeed  and  c is  the  speed  of  sound) . The 

3 3 3 

time  t is  referenced  to  the  time  of  the  midpoint  of  the  0.5-s  sample  of  data 
R 

that  was  digitized  to  determine  the  set  of  measured  input  sound  pressure 
levels. 

If  the  time  at  overhead  (TOH  at  <f  * 90°)  and  the  time  at  the  start  of 
digitizing  the  0.5-s  samples  of  data  (TS)  are  known,  and  the  time  of  the 
selected  noise  data  sample  relative  to  the  start  time  (TN)  is  also  known, 
then  the  relative  time  tj^  can  be  found  from 

t_  - TS  + 0.25  + TN  - TOH  (17) 

K 

in  seconds  if  TS  and  TOH  are  in  seconds  on  a consistent  basis.  The  0.25-s 
term  is  Included  to  shift  the  origin  to  the  middle  of  the  0.5-s  period  since 
TS  starts  at  the  beginning  of  the  first  data  sample  for  this  demonstration. 
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From  geometrical  considerations,  the  time  tj^  can  be  shown  to  be  related 
to  the  height,  airspeed,  Mach  number,  and  propagation  angle  through 

t = (AMH/V  ) (M  CSC  ip  - cot  ((»)  • (18) 

K d d 

Note  that  for  t„  = 0,  equation  (18)  requires  that  cos  = M . 

K 3 

We  solve  equation  (18)  for  cos  ip  using  some  algebraic  manipulation  and 
trigonometric  Identities  to  find 

cos  = [M  (AMH/V  t_)^]  [(AMH/V  t„)^  + 1]~^ 

3 3 K 3 K 

+ [(AMH/V^tj^)^  + 1]“^ 

jMf(AMH/V^t_)^  - [(AMH/V  t_)^  + 1]  [(AMH/V  t_)V  - 1]'  (19) 

^3  3rv  3K  31\3 

where  the  + sign  applies  when  t_  < 0 (i.e.,  for  0®  < < cos  ^ M ) and  the  - 

-1  ^ 

sign  applies  when  tj^  > 0 (I.e.,  for  cos  M^  < ij;  < 180°). 

With  a value  for  ip  from  equation  (19),  the  total  length  of  the  propaga- 
tion path  is 

PD  - AMH/sln  ip.  (20) 

The  length  D of  any  segment  along  the  propagation  path  defined  by  heights 

J I 

HM  where  meteorological  data  were  measured  can  be  determined.  In  general,  from  | 

Dj  - - HMj)/sin  ip  (21) 

except  for  the  first  and  last  segments. 

For  the  first  segment,  we  sort  through  the  HM  array  to  find  the  first 
height  value  that  is  greater  than  the  microphone  height  ZM  and  call  that  value 
HMjj.  Then,  the  length  of  the  first  segment  is 

- (HMjj  - ZM)/sin  ip.  (22) 

For  the  last  segment,  we  sort  through  the  HM  array  and  find  that  height 
value  which  is  just  less  than  the  airplane  height  AH  and  call  that  value 
The  length  of  the  last  segment  is  then 

- (AH  - HMj^)/sin  4».  (23) 


24 


The  Input  data  described  in  figure  3 are  read  or  calculated.  With  the 
array  of  propagation  distances  computed  using  equations  (21),  (22),  and  (23), 
calculation  of  atmospheric  absorption  adjustments  can  proceed. 

Two  test  cases  were  selected.  The  measured,  or  test,  meteorological 
data  for  the  two  cases  are  shown  in  figure  5 as  vertical  profiles  of  air 
temperature  and  relative  humidity.  There  were  31  measurements  at  heights 
between  1.2  and  915  m.  The  station  barometric  pressure  BP  is  given  in  the 
subcaptions.  Reference  profiles  of  temperature  and  relative  humidity  are 
shown  by  the  dashed  lines.  The  airplane  height  was  629  m for  test  case  1 
and  154  m for  test  case  2 as  indicated  by  the  arrows  marked  AH.  Winds  were 
always  calm. 

For  test  case  1 the  temperature  above  200  m was  not  a great  deal 
different  than  the  reference  temperature,  below  200  m the  air  was  colder 
than  the  reference  temperature.  Relative  humidity  was  everywhere  significantly 
less  than  the  reference  relative  humidity. 

For  test  case  2,  the  temperature  was  everywhere  less  than  the  reference 
temperature,  especially  below  200  m.  The  relative  humidity  was  close  to  the 
reference  value  above  150  m and  greater  than  the  reference  value  below  150  m. 

For  test  case  1,  the  molar  concentration  of  water  vapor  ranged  from  1.0 
to  0.45  percent  as  height  increased  and  from  1.03  to  0.93  percent  for  test 
case  2.  Thus  the  air  was  always  significantly  drier  than  at  the  reference 
condition  of  25°C  and  70  percent  relative  humidity  at  10  m where  the  molar 
concentration  of  water  vapor  is  approximately  2.2  percent.  Atmospheric 
absorption  under  the  test  conditions  would  therefore  be  expected  to  be  sub- 
stantial for  the  higher  frequencies  and  the  adjustments  to  reference  condi- 
tions would  be  large  for  the  higher  frequency  bands,  especially  over  the  long 
propagation  distance  of  test  case  1. 

The  results  of  running  program  BANDLOSS  for  test  case  1 are  shown  in 
table  1 in  the  format  provided  by  the  program.  Table  2 lists  the  temperature 
and  relative  humidity  profile  data  for  test  case  1.  Tables  3 and  4 provide 
similar  data  for  test  case  2. 
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1 Table  1.-Results  from  program  BANDLOSS  for  test  case  1. 
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AH  = airplane  height  above  ground  level;  AS  = airspeed;  AM  = airplane  Mach  number;  ZM  = height  of 
microphone  above  ground  level;  PD  = sound  propagation  distance;  TS  = time  for  start  of  data  analysis; 
TN  ° time  after  TS  for  noise  data  sample;  TOH  • time  when  airplane  was  over  microphone;  and  PSI  = 


sound  propagation  angle  relative  to  flight  path. 
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Table  2.-Measured  and  reference  profile  data  for  air  teniperature 
and  relative  humidity  applicable  to  test  case  1. 
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® 'Measured'  temperatures  and  relative  humidities  listed  above  are  values  interpolated  from  the  actual  measured 
values  to  represent  temperatures  and  relative  humidities  at  heights  midway  between  the  heights  where  the 
temperatures  and  relative  humidities  were  actually  measured.  The  midpoint  values  represent  average  values 
over  the  pathlength  segments. 
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Table  3.-Resuit«  from  program  BANDLOSS  for  test  cate  2. 
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Table  A.-Measured  and  reference  profile  data  for  air  temperature 
and  relative  humidity  applicable  to  text  cate  2. 
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The  results  for  the  test  cases  are  shown  graphically  In  figure  6.  The  j 

subcaptions  give  the  airplane  height,  airspeed,  Mach  number,  microphone 
height,  time  of  the  midpoint  of  the  data  sample  relative  to  the  overhead  time, 
total  propagation  distance,  and  sound  emission  angle.  Test  case  1 was  for  a j 

time  prior  to  overhead;  test  case  2 was  for  a time  after  overhead. 

For  test  case  2,  all  24  band  levels  represent  actual  measured  data.  For  , 

test  case  1,  however,  the  levels  in  the  last  five  bands  were  provided  by  i 

extrapolation  and  thus  the  input  spectrum  for  test  case  1 is  partially  l 

artificial.  ■ 

1 

Referring  to  tables  1 and  3,  we  note  that  the  calculated  adjustment  j 

factors  are  actually  negative  numbers  for  several  of  the  'midfrequency'  bands  | 

—630  to  1250-Hz  for  test  case  1 and  400  to  3150-Hz  for  test  case  2.  The  i 

negative  values  indicate  less  atmospheric  absorption  under  the  test  conditions  ■ 

than  under  the  reference  conditions,  in  contrast  to  the  usual  situation  where  j 

i 

the  minimum  absorption  occurs  near  the  reference  conditions. 

In  the  low  frequency  bands,  the  adjustment  factors  are  all  small  to 
negligible,  as  expected,  even  for  the  899-m  pathlength  of  test  case  1. 

In  the  high-frequency  bands,  the  adjustment  factors  are  quite  large, 
especially  for  test  case  1.  For  test  case  2,  the  adjusted  spectrum  shape 
looks  plausible.  For  test  case  1,  however,  the  relatively  flat  spectrum  of 
the  adjusted  sound  pressure  levels  may  be  unreasonable,  although  its  appearance 
is  distorted  in  figure  6(a)  because  of  the  need  to  use  a 20  dB/division  scale 
for  test  case  1 instead  of  the  standard  10  dB/dlvision  scale  as  used  for 
test  case  2 in  figure  6(b). 

For  test  case  1,  moreover,  the  last  five  band  levels  (4000  to  10,000  Hz) 
had  to  be  estimated  as  noted  earlier  because  the  measurement  system  noise 
floor  was  equivalent  to  a sound  pressure  level  of  about  40  dB.  The  true 
spectral  rolloff  rate  might  very  well  have  been  more  rapid  than  estimated  for 
the  meteorological  conditions  of  figure  5(a)  and  the  propagation  distance  of 
899  m.  If  a more  rapid  rolloff  had  been  estimated,  the  calculated  values  of 
the  adjusted  sound  pressure  levels  would  have  been  lower.  More-rapid  rolloff 
rates,  however,  would  probably  have  been  no  more  accurate  than  the  arbitrary 
rates  selected  for  test  case  1. 
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(a)  TEST  CASE  1;  AH  = 629.0  m,  AS  = 81.5  m/s,  AM  = 0.24, 
ZM  = 1.2  m,  TR  = -5.25  s,  PD  = 899.1  m,  PSI  = 44.3° 


1/3-OCTAVE-BAND  CENTER  FREQUENCY,  kHz 

(b)  TEST  CASE  2:  AH  = 154.0  m,  AS  = 74.4  m/s,  AM  = 0.22, 

ZM  = 1.2  m,  TR  = +5.75  s,  PD  = 377.2  m,  PSI  = 156.1° 

Figure  6.-ResutU  of  using  program  BANDLOSS  to  determine  atmospheric  absorption 
adjustments  from  test  to  reference  meteorological  conditions. 
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Test  case  1 was  included  deliberately  even  though  some  band  levels  were 
missing  and  had  to  be  supplied  by  arbitrary  extrapolation.  The  results  indi- 
cate the  difficulty  of  obtaining  reasonable  estimates  of  the  true  pressure 
' , spectral  density  at  the  microphone  on  the  basis  of  band  levels  when  the 

: spectral  slopes  are  large. 

Large  spectral  slopes  likely  mean  that  the  single  straight-line  approxima- 
tion of  Appendix  B is  not  realistic  for  the  pressure  spectral  density  function 
over  the  entire  passband  of  the  filters.  Large  spectral  slopes  probably  also 
[ mean  that  significant  energy  is  contributed  by  either  the  lower  or  upper 

stopbands.  Thus  the  fundamental  assumption  that  the  actual  frequency  response 
I of  the  real  filters  can  be  reasonably  approximated  by  that  of  ideal  filters  is 

(■  probably  not  valid  for  spectral  slopes  as  steep  as  those  for  the  estimated 

i 

i high-frequency  sound  pressure  levels  of  test  case  1.  If  the  filter  response 

cannot  be  represented  by  that  of  ideal  filters,  then  the  analysis  of  Appendix 
B would  have  to  be  modified  to  include  integration  over  that  portion  of  the 
filter's  stopbands  which  contribute  'significant'  energy.  A definition  of 
'significant'  would  have  to  be  supplied  and  a test  for  convergence  of  the 
integral  would  have  to  be  developed. 

The  conclusions  are  (1)  that  extrapolation  to  supply  missing  low-  or 
high-frequency  band  levels  may  give  misleading  results,  especially  extrapola- 
tion to  obtain  high-frequency  sound  pressure  levels  where  atmospheric  absorp- 
tion effects  are  largest,  and  (2)  that  there  probably  is  some  maximum  spectral 
slope  for  which  the  assumptions  of  Appendix  B are  valid.  It  was  beyond  the 
scope  of  the  project  to  determine  the  range  of  spectral  slopes  for  which  the 
method  of  Appendix  B was  valid. 

DATA  PROCESSING  AND  ANALYSIS 

Most  existing  systems  for  processing  a high-volume  of  aircraft  flyover 
noise  measurements  should  already  have  provisions  to  supply  almost  all  the 
input  data  needed  by  program  BANDLOSS  for  calculating  atmospheric  absorption 
loss  adjustments  for  bands  of  noise.  Some  minor  modifications  may  be  needed 
to  provide  data  in  the  measurement  units  required  by  the  program.  A routine 
may  be  required  to  calculate  the  lengths  of  the  sound  propagation  paths  as  is 


differences  (DL  in  lines  176  to  179)  and  the  definitions  for  RF,  SE,  and 
ISBN  would  have  to  be  added.  The  calculations  of  reference  meteorological 
conditions  would  also  have  to  be  added. 

The  statements  in  lines  191  to  233,  with  subroutine  ANSAB  (lines  278  to 

297),  comprise  the  essence  of  the  program  statements  needed  to  determine  band- 

loss  adjustment  factors  AL„  The  results  of  the  calculations  could  be 

D I 1 K 

applied  to  any  0.5-s  sample  of  wide-band  data  included  in  the  analysis  of  a 
flyover  noise  signal. 

The  three  special  considerations  listed  earlier  should  be  taken  into 
account  before  attempting  to  calculate  absorption  adjustments  to  measured 
sound  pressure  levels.  These  were:  (1)  determining  which  of  the  input  sound 
pressure  levels  are  missing  because  of  interference  from  high  ambient  or 
background  noise  levels,  (2)  removing  spectral  irregularities  caused  by 
ground  reflection  effects,  and  (3)  separating  the  input  spectrum  into  discrete 
and  broadband  components. 

Special  rules  will  have  to  be  developed  to  determine  reasonable  estimates 
for  the  slope  of  the  input  spectrum  when  it  is  determined  that  there  are 
missing  band  levels  between  the  first  and  last  bands.  A separate  computational 
routine  would  have  to  be  created  to  compute  the  adjustment  factors  for  the 
discrete-frequency  components  of  the  spectrum,  AL^  based  on  the  procedures 
of  Appendix  A. 

RECOMMENDATIONS 

This  report  has  described  the  development  of  a method  to  calculate  the 
atmospheric  absorption  experienced  by  broadband  sound  propagating  in  a 
quiescent  atmosphere  and  analyzed  by  fractional  octave  band  filters.  The 
analytical  method  for  bands  of  noise  is  based  on  several  assumptions.  Some 
additional  analyses  and  experimental  data  are  needed  to  verify  the  assump- 
tions before  adopting  the  procedures  in  routine  applications. 

The  recommendations  given  below  are  divided  into  three  categories: 
general,  procedural  and  instrumentation  considerations,  and  validation. 
Recommendations  that  are  general  in  nature  are  for  studies  to  help  develop 
the  method  for  use  in  a variety  of  applications,  such  as  research  investiga- 
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tions  of  aircraft  noise  propagation  over  long  distances  or  through  atmos- 
pheric turbulence,  prediction  of  aircraft  flyover  noise  levels  or  noise  ex- 
posure level  contours,  and  aircraft  noise  type  certification.  The  procedural 
and  instrumentation  considerations  are  relevant  to  all  applications,  but  are 
specifically  appropriate  to  the  aircraft  noise  measurement  procedures  of 
Appendix  A of  FAR  36. 

General 

1.  One  of  the  fundamental  assumptions  in  the  analysis  is  that  the  slope 
of  the  pressure  spectral  density  y£  frequency  curve  for  the  sound  incident  on 
a microphone  is  not  'too'  steep.  This  assumption  permits  replacement  of  the 
frequency  response  of  the  actual  filters  used  in  processing  the  signal  by  the 
frequency  response  of  ideal  filters. 

A study  should  be  conducted  to  determine  the  range  of  positive  and  nega- 
tive spectral  slopes  for  which  the  assumption  of  ideal  filters  gives  valid 
results.  The  study  should  consider  the  response  characteristics  of  various 
actual  analog  and  digital  filters  as  well  as  the  response  characteristics  of 
filters  meeting  the  minimum  transmission  loss  characteristics  defined  in 
reference  12  for  octave  and  1/ 3-octave-band  filters. 

2.  The  magnitude  of  the  calculated  absorption  depends  on  the  spectrum 
of  the  sound  incident  on  the  microphone.  The  actual  sound  spectrum  is 
approximated  by  a series  of  straight-line  segments.  A fundamental  assumption 
is  that  the  slopes  of  the  straight-line  segments  can  be  determined  from  the 
average  sound  pressure  spectral  densities  at  the  center  frequencies  of  the 
filter  bands.  For  ideal  filters,  the  bandwidth  of  the  filters  is  defined 
precisely.  For  real  filters,  an  appropriate  definition  of  bandwidth  must 

be  developed. 

The  analysis  in  Appendix  B develops  a procedure  for  estimating  the  slopes 
of  the  straight-line  segments  approximating  the  true  spectrum  on  the  basis  of 
the  difference  in  the  average  pressure  spectral  densities  in  the  frequency 
bands  above  and  below  the  particular  band  of  interest.  The  validity  of  the 
procedure  needs  to  be  critically  examined,  especially  for  measured  sound 
spectra  which  have  large  sound  pressure  level  differences  in  adjacent  bands 
as  often  encountered  in  high-frequency  data.  The  examination  should  be  con- 
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ducted  In  conjunction  with  the  study  of  the  validity  of  approximating  the 
response  of  real  filters  by  that  of  ideal  filters.  The  impact  of  alternate 
schemes  for  estimating  the  spectral  slopes  should  be  evaluated  for  typical 
spectra  produced  by  a variety  of  aircraft  types. 

3.  With  the  assumption  of  ideal  filters,  the  procedure  for  determining 
the  atmospheric  absorption  loss  of  broadband  noise  includes  evaluation  of  an 
integral  whose  limits  are  the  lower  and  upper  bandedge  frequencies  of  the 
corresponding  ideal  fractional-octave-band  filter.  The  integration  is  per- 
formed by  dividing  the  fractional  octave  band  into  several  smaller  frequency 
bands  and  approximating  the  integrand  by  a straight  line  (on  logarithmic 
scales)  over  each  sub-band.  The  examples  in  this  report  considered  1/3-octave 
bands  which  were  then  each  subdivided  into  twenty  1/60-octave  bands. 

The  accuracy  of  the  Integration  Improves  as  the  fractional  octave  band  is 
subdivided  into  an  increasing  number  of  sub-bands.  The  computing  time, 
however,  also  increases.  A study  should  be  performed  to  evaluate  the  sensiti- 
vity and  accuracy  of  the  results  to  the  number  of  sub-bands  in  the  integra- 
tion for  a variety  of  assumed  meteorological  conditions  and  sound  spectra. 

4.  In  the  demonstrations  of  the  applications  of  the  procedure  developed 
here,  a key  assumption  was  that  the  noise  produced  by  an  aircraft  flying 
over  a microphone  near  the  ground  originated  from  a single  acoustic  point 
source.  The  source  was  assumed  to  be  located  at  a point  on  the  airplane 
flight  path  determined  by  the  time  of  emission.  The  assumption  of  an 
acoustic  point  source  is  equivalent  to  assuming  that  the  microphone  is  in  the 
geometric  and  acoustic  far  field  at  all  frequencies  of  interest. 

For  many  cases  of  practical  importance,  the  assumption  of  a single 
equivalent  point  source  may  not  be  reasonable  because  the  airplane  is  a large, 
multi-engine  type.  The  minimum  distance  between  the  microphone  and  the  air- 
plane reference  point  may  be  relatively  short.  Specific  procedures  need  to  be 
developed  for  application  in  situations  where  the  aircraft  cannot  be  consider- 
ed to  be  equivalent  to  a single  acoustic  point  source.  The  procedures  should 
consider  the  necessity  of  defining  apparent  sound  source  locations  as  a 
function  of  frequency.  For  example,  low-frequency  jet  noise  sources  might 
be  located  farther  behind  the  jet  exhaust  nozzle  than  high-frequency  jet  noise 
sources. 
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5.  The  sound  propagation  path  was  assumed  to  be  a straight  line  between  j 

the  location  of  the  effective  acoustic  point  source  and  the  microphone. 

Bending  or  curvature  of  the  sound  propagation  path  because  of  gradients  In 
temperature  or  wind  was  ignored  even  though  gradients  were  included  In  the 
profiles  of  test  and  reference  temperature.  A study  should  be  performed  to 
evaluate  the  need  to  incorporate  a ray-tracing  procedure  in  the  determination 
of  the  lengths  and  directions  of  the  various  segments  of  the  sound  propaga- 
tion path. 

6.  The  demonstrations  in  this  report  of  the  procedure  for  calculating 
atmospheric  absorption  assumed  that  the  nonlinear  effects  associated  with 
propagation  of  high-amplitude  sound  waves  could  be  neglected.  It  is  known, 
however,  that  the  spectrum  of  high-amplitude  sound  waves  is  modified  by  non- 
linear effects  as  the  waves  propagate  away  from  the  source. 

The  sources  of  aircraft  noise  are  obviously  sources  of  high-amplitude 
sound.  Research  is  needed  to  quantify  the  spectral  effects  introduced  by 
nonlinear  propagation  of  high-amplitude  sound  waves.  Procedures  should  then 
be  developed  to  incorporate  these  spectral  effects  into  the  procedures  for 
determining  atmospheric  absorption  adjustments.  Since  nonlinear  effects 
modify  the  shape  of  the  sound  spectrum  and  since  spectral  slope  is  a factor 
in  determining  the  atmospheric  absorption  of  bands  of  noise,  some  method  of 
simultaneously  accounting  for  both  effects  may  be  required. 

Procedural  or  Instrumentation  Considerations 

1.  To  determine  an  atmospheric  absorption  factor  to  adjust  measured 
sound  pressure  levels  from  test  to  reference  meteorological  conditions,  it 
was  necessary  to  define  an  acoustic  reference  atmosphere  with  reference  lapse 
rates  for  temperature,  relative  humidity,  and  pressure.  Reference  meteorolo-  ! 

gical  conditions  at  a height  of  10  m above  mean  sea  level  (altitude)  were  J 

specified  to  be  those  defined  in  FAR  36. 

It  is  recommended  that  the  acoustic  reference  atmosphere  defined  in  this 
report  be  adopted  for  the  purpose  of  calculating  atmospheric  absorption  ad- 
justment factors. 
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2.  Since  the  atmospheric  absorption  correction  for  bands  of  noise  depends 
on  the  sound  spectrum  which  must  be  approximated  from  the  measured  sound 
pressure  levels,  a problem  arises  when  some  of  the  measured  band  levels  are 
missing  as  a result  of  interference  from  high-level  background  noise. 

Special  rules  for  handling  the  problem  of  missing  band  levels  need  to 
be  developed  and  tested.  The  rules  should  provide  a procedure  for  examining 
the  set  of  measured  sound  pressure  levels  and  either  estimating  an  appropriate 
spectral  slope  or  deciding  that  no  rational  estimate  can  be  made  for  the 
slope  and  that  it  is  not  feasible  to  determine  corrections  for  differences  in 
atmospheric  absorption  for  that  particular  data  sample. 

3.  A procedure  to  examine  the  measured  set  of  sound  pressure  levels  and 
remove  ground-reflection  spectral-interference  effects  should  be  developed, 
evaluated,  and  incorporated  with  the  atmospheric  absorption  procedure.  The 
evaluation  should  determine  the  magnitude  of  differences  between  removing  or 
not  removing  the  ground-reflection  effects  before  calculating  adjustments  for 
atmospheric  absorption.  Microphone  heights  of  1.2  and  10  m should  be  evaluat- 
ed. The  analytical  procedure  for  removing  ground  reflection  effects  could  be 
based  on  the  techniques  of  references  13  and  14. 

4.  Since  different  procedures  are  used  to  calculate  atmospheric  absorp- 
tion losses  for  the  discrete-frequency  and  the  broadband  components  of  an 
aircraft  noise  signal,  a procedure  for  separating  a given  set  of  measured  air- 
craft noise  levels  into  the  discrete  and  broadband  components  needs  to  be  in- 
corporated into  an  overall  procedure  for  determining  atmospheric  absorption 
adjustments.  The  pure-tone  identification  procedure  of  FAR  36  would  be  a 
logical  candidate  for  performing  the  spectral  separation.  Samples  of  aircraft 
noise  signals  with  a mixture  of  discrete  and  broadband  components  should  be 
evaluated  to  test  the  combined  procedure. 

5.  The  specific  time  (e.g.,  beginning,  midpoint,  or  end)  for  synchroniz- 
ing the  0.5-s  data  sample  with  the  aircraft's  position  is  a factor  in  determin- 
ing the  sound  propagation  distance.  Data  averaging  time  and  the  data  averag- 
ing procedure  are  factors  to  be  considered  in  relating  the  measured  sound 
pressure  levels  to  the  average  absorption  experienced  by  sound  waves  incident 
on  a microphone  over  the  duration  of  a given  data  sample.  For  consistent 
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determination  of  atmospheric  absorption  adjustments,  it  is  recommended  that  the 
midpoint  of  the  duration  of  the  data  sample  be  used  for  synchronizing  acoustic 
and  aircraft  data  and  that  linear  integration  over  the  data-sample  duration  be 
used  to  obtain  the  time  average  of  the  filtered  signals. 

6.  The  current  version  of  FAR  36  requires  application  of  atmospheric 
absorption  adjustments  to  only  one  of  the  0.5-s  data  samples,  the  one 
associated  with  PNLTM^gg^,  in  order  to  determine  PNLTMj.g£.  For  most  research 
and  prediction  purposes,  it  is  recommended  that  absorption  adjustments  be 
determined  for  each  data  sample  throughout  the  total  duration  of  Interest. 

For  aircraft  noise  type  certification,  it  is  also  recommended  that  considera- 
tion be  given  to  requiring  calculation  of  atmospheric  absorption  adjustments 
for  all  data  samples  within  the  duration  of  Interest. 

A study  should  be  conducted  to  evaluate  the  Impact  of  changing  from 
applying  absorption  adjustments  to  just  one  data  sample  to  applying  them  to 
all  relevant  data  samples.  The  study  should  consider  the  impact  on  aircraft 
noise  signals  from  a variety  of  aircraft  types  at  various  engine  power  settings 
and  distances  and  under  various  test  meteorological  conditions.  One  purpose 
of  the  study  would  be  to  determine  if  the  recommended  procedure  shifts  the 
time  of  occurrence  of  PNLTM  under  the  acoustic  reference  meteorological  con- 
ditions to  a different  time  than  under  the  test  meteorological  conditions. 

It  is  also  recommended  that  the  study  evaluate  the  impact  of  determining 
the  duration  correction  factor  under  acoustic  reference  meteorological  con- 
ditions from  the  set  of  PNLTj-gf  values,  calculated  from  the  series  of  adjusted 
sound  pressure  levels,  at  the  times  corresponding  to  the  PNLTj-gf  values  which 
are  10  dB  less  than  the  value  of  PNLTMj.g£.  The  current  FAR  36  procedure  of 
accounting  for  atmospheric  absorption  for  just  one  data  sample  does  not 
specifically  consider  the  change  in  duration  that  would  be  associated  with  the 
time  variation  of  tone-corrected  perceived  noise  level  under  acoustic  reference 
meteorological  conditions. 

Validation 

1.  The  complexity  and  Importance  of  the  problem  of  adjusting  measured 
data  from  test  to  reference  conditions  require  that  the  total  impact  of  the 
proposed  new  procedures  be  carefully  evaluated.  The  evaluation  should  consider 
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the  Impact  on  the  noise  produced  by  different  types  of  jet  and  propeller-powered 
aircraft  tested  at  various  power  settings  and  distances  and  under  a variety  of 
meteorological  conditions.  Tests  under  meteorological  conditions  which  are 
close  to  those  of  the  acoustic  reference  atmosphere  should  be  conducted  to 
provide  the  baseline  data  against  which  to  compare  the  adjusted  data  from  tests 
under  other  meteorological  conditions. 

Differences  between  test-day  and  acoustic-reference-day  sound  pressure 
levels,  perceived  noise  levels  (with  and  without  tone  corrections),  and 
effective  perceived  noise  levels  should  be  evaluated.  Atmospheric  absorption 
adjustment  procedures  should  include  the  method  described  in  this  report  and 
the  method  in  FAR  36  which  is  based  on  that  in  SAE  ARP  866A. 

2.  As  part  of  the  overall  validation  effort,  it  is  recommended  that  ^ 

consideration  be  given  to  determining  whether  the  proposed  new  procedure  can 
be  .simplified  without  significant  loss  of  accuracy. 

With  the  procedure  of  ARP  866A,  there  are  essentially  only  four  variables: 
distance,  band  center  frequency  (or  lower  bandedge  frequency  for  the  high- 
frequency  bands),  air  temperature,  and  relative  humidity.  For  pure  tones,  the 
proposed  new  procedure  has  these  same  four  variables  (but  with  the  measured  j 

tone  frequency  rather  than  band  center  frequency) , plus  the  pressure  of  the 
air.  j 

For  bands  of  noise,  the  proposed  new  procedure  requires  integration  over 
the  range  of  frequencies  contributing  significant  energy  to  the  indicated 

band  level.  The  integration  process  requires  determination  of  a rational  j 

estimate  of  the  slope  of  the  sound  spectrum  over  the  same  range  of  frequency. 

If  the  assumption  of  ideal  filters  does  not  turn  out  to  be  acceptable,  then 
the  integration  may  also  require  incorporation  of  the  response  of  the  actual 
filters  used  for  analyzing  the  measured  sound  pressure  levels. 

With  only  the  variables  of  frequency,  temperature,  and  relative  humidity, 
it  was  feasible  for  SAE  ARP  866A  to  Include  tables  showing  the  dependence  of 
sound  absorption  coefficients  on  the  variables.  With  the  proposed  new  pro- 
cedure, it  is  not  feasible  to  construct  tables  for  the  sound  absorption 
coefficient  unless  the  ranges  of  some  of  the  variables  are  greatly  restricted. 
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One  potential  simplication,  however,  that  should  be  Included  in  the  validation 
effort  is  replacement  of  the  integration  by  determination  of  the  differences 
in  atmospheric  absorption  at  the  band  center  frequency  only, 

' CONCLUDING  REMARKS 

, 1.  A computational  procedure  and  corresponding  computer  program  were 

developed  to  determine  differences  in  atmospheric  absorption  under  test  and 

I reference  meteorological  conditions.  The  procedure  is  based  on  calculation  of 

pure-tone  atmospheric  absorption  according  to  the  experimentally  verified 

I method  in  the  proposed  American  National  Standard  ANS  SI. 26.  The  absorption 

adjustment  procedure  developed  here  is  applicable  to  discrete-frequency  and 

I 

i broadband  components  of  an  aircraft  noise  signal. 

2.  An  acoustic  reference  atmosphere  was  defined  for  the  variation  of 
temperature,  relative  humidity,  and  pressure  with  height  above  mean  sea 
level.  Reference  conditions  from  FAR  36  were  specified  at  a height  of  10  m. 

3.  Application  of  the  data  adjustment  procedure  was  demonstrated  for  two 
sets  of  1/ 3-octave-band  aircraft  noise  levels.  The  adjustment  procedure  for 
bands  of  noise  requires  evaluation  of  an  integral  over  the  frequency  range  of 
each  band.  To  perform  the  integration,  an  estimate  must  be  supplied  of  the 
slope  of  the  spectrum  of  the  sound  pressure  Incident  on  the  microphone.  The 
computational  procedure  includes  a method  for  estimating  the  spectrum  slope 
on  the  basis  of  the  measured  1/ 3-octave-band  sound  pressure  levels. 

4.  Recommendations  were  developed  for  several  general  and  specific  studies 
to  provide  information  for  validating  or  refining  the  assumptions  made  in 
developing  the  analytical  procedure  for  calculating  absorption  adjustments  for 
bands  of  noise.  Other  studies  were  recommended  as  part  of  a validation  effort 
that  should  be  conducted  prior  to  considering  Incorporation  of  the  proposed 
new  procedure  in  FAR  36  as  a replacement  for  the  existing  procedure  that  is 
based  on  SAE  Aerospace  Recommended  Practice  ARP  866A. 

DyTec  Engineering,  Inc. 

Huntington  Beach,  California  92649 
30  December  1977 
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ATMOSPHERIC  ABSORPTION  ADJUSTMENT  FOR  PURE  TONES 

In  this  Appendix  we  develop  the  analytical  expressions  needed  to  calcu- 
late the  magnitude  of  the  atmospheric-absorption  adjustment  factor  between 
test  and  reference  conditions  for  discrete-frequency  pure-tone  signals. 
Absorption  losses  are  determined  in  accordance  with  the  procedures  of  the  pro- 
posed American  National  Standard  ANS  Si. 26. 

The  amplitude  of  the  pressure  in  a sound  wave  decreases  exponentially  in 
the  direction  of  propagation  because  of  absorption  processes.  References  1 
and  2 describe  the  physical  mechanisms  for  these  processes. 

Ignoring  attenuation  of  a sound  wave's  amplitude  caused  by  spreading 
losses  as  the  distance  from  the  source  increases  and  assuming  a quiescent 
atmosphere,  the  decrease  in  pressure  from  point  1 to  point  2 over  a distance 
C can  be  expressed  as 

P2  = Pi  e (Al) 

where  a is  the  sound  attenuation  coefficient  due  to  atmospheric  absorption, 
or,  for  short,  the  absorption  coefficient  in  units  of  nepers/m  when  distance 
C is  in  meters. 

The  absorption  coefficient  a is  a function  of  the  frequency  f of  the 
sound  wave  and  the  temperature,  humidity,  and  pressure  of  the  air.  From  the 
proposed  ANS  SI. 26,  the  formula  giving  the  relation  between  a and  the  physical 
variables  is 

a = f^  {1.84-10"^^  (T/Tq)^''^ 

+ (T/TQ)“^^^tl.278-10‘^(exp(-2239.1/T)/(f^^Q  + 

+ 0.1068(exp(-3352/T))/(f^  ^ + (f^/f^  j,))]}  (A2) 

where  P is  atmospheric  pressure,  Pq  a reference  atmospheric  pressure,  T the 
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temperature  of  the  air,  Tq  a reference  air  temperature,  f^  ^ the  vibrational 
relaxation  frequency  for  oxygen,  and  f the  vibrational  relaxation  fre- 
quency  for  nitrogen. 

The  expressions  for  the  two  vibrational  relaxation  frequencies  are 

fr  0 = (P/Po)^24  + 4.41-10\  [(0.05  + h)/(0.391  + h)]}  (A3) 

and 

^r  N “ (P/Po)(T/Tq)"^^^[9  + 350h  exp{-6. 142 [ (T/Tq)“^^^  - 1]}]  (A4) 

where  h is  the  molar  concentration  of  water  vapor  in  a given  sample  of  moist 
air. 

The  value  of  the  molar  concentration  is  determined  from  the  relative 
humidity  and  atmospheric  pressure  using 

" ■ "t 

where  h^  is  the  relative  humidity  and  P^^^  is  the  saturation  vapor  pressure 
of  water  (over  liquid  water)  at  the  pressure  and  temperature  of  the  moist  air. 

The  ratio  P /P_  is  a function  of  temperature  only  and,  as  in  proposed 
sat  u 

ANS  SI. 26,  is  determined  from 

log  (Pg^t^Po)  “ 10.79586[1  - (Tqj^/T)] 

- 5.02808  log  (T/Tqj^) 

+ 1.50474.l0-\l  - io-«- 29^92 [(T/T^jP  -1]^ 

+ 0.42873*10"^(10''*'^^^^^^^  ■ ('^01^^^^  - 1) 

-2.2195983  (A6) 

where  T^^  is  the  triple-point  isotherm  temperature  with  the  exact  value  of 
273.16  kelvlns. 

Equations  (A2)  to  (A6)  are  combined  to  determine  a value  for  a for  given 
values  of  f,  T,  h^  and  P and  suitable  values  for  Pq  and  T^.  In  the  procedure 
described  in  this  report,  values  of  a are  calculated  by  a subroutine  called 
ANSAB  for  American  National  Standard  ABsorptlon.  Reference  pressure  and 
temperature  values  are  Pq  = 101.325  kPa  and  Tq  =■  293.15  K. 
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Over  any  given  path  of  length  C where  the  absorption  has  the  value  a,  the 
change,  in  decibels,  in  pure-tone  sound  pressure  level  AL^  between  that  at  the 
beginning  and  end  of  the  path  (in  the  direction  of  propagation)  is  determined 
from  the  squared  pressures  and  equation  (Al)  as 

r 

AL^  = - -10  log  (P2/Pj)  = -10  log  e"^“^  (A7) 

= 20  aC  log  e 

= 8.686aC  (A8) 

= aC 

where  a = 8.68b  a.  is  the  pure-tone  absorption  coefficient  in  units  of  decibels 
per  meter  or  dB/m  and  the  symbol  log  stands  for  logarithm  to  the  base  10. 

[Note  that  a would  have  units  of  dB/(100m)  if  C were  in  hundreds  of  meters 
or  units  of  dB/km  if  f were  in  kilometers.]  Subroutine  ANSAB  also  determines 
values  for  the  absorption  coefficient  a. 

If  the  temperature,  relative  humidity,  and  pressure  are  not  constant  over 
the  propagation  distance  C,  then  the  pure-tone  absorption  loss  should  be 


calculated  from 


AL^  = 8.686  J ^ a di 
T 


(AlO) 


where  and  C2  bhe  coordinates  at  the  beginning  and  end  the  propaga- 
tion path. 

With  AC  = Co  - Cl  and  an  average  absorption  coefficient  a defined  by 
— C ^ 

a = (1/AC)  £-/  ^ n dC,  equation  (AlO)  can  be  written  as 


AL.J,  *•  8.686  a AC 


where  a is  the  mean  value  of  a over  pathlength  AC. 


(All) 


To  determine  the  difference  in  atmospheric  absorption  loss  between  test 
and  reference  meteorological  conditions  along  a sound  propagation  path,  we  let 
in  equation  (A7)  represent  the  sound  pressure  level  at  the  receiver  after 
passage  through  the  test  atmosphere,  or  atmosphere  1.  The  level  L2  is  the 
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sound  pressure  level  at  the  receiver  (from  the  same  source  at  the  same 
location)  after  passage  through  the  reference  atmosphere,  or  atmosphere  2. 

The  pure-tone  atmospheric-absorption  adjustment  factor  between  test  and 

reference  conditions,  AL  , Is,  by  analogy  to  equation  (A7), 

1 f IK 


^^T,TR  ^2  H ^Pr2'^Pr1^ 


(A12) 


2 2 

where  p^^  and  Pj^2  the  squared  pure-tone  sound  pressures  at  the  receiver 
and  AL„  _ Is  the  number  of  decibels  to  be  added  algebraically  to  the  measured 

1 y 1 K 


levels,  to  get  the  adjusted  levels  L2. 


If  the  squared  sound  pressure  at  the  source  Is  Pg(f,0)  and  If  the  propa- 
gation path  consists  of  segments  A^  long,  then  the  squared  sound  pressure 
2 

at  the  receiver,  p^, , can  be  written  by  repeated  application  of  equation  (Al), 
RJ. 

as 

-2a  AC,  -2a  AC,  -2a  AC„ 

PrI  ° [Pg(f,0)](e  ^)(e  •••  (e  ")  (A13) 

where  the  subscripts  Indicate  the  path  segment  number  and  the  conditions  of 
atmosphere  1.  There  are  n segments  along  the  total  length  of  the  path.  The 
absorption  coefficients  would  be  those  calculated  using  equations  (A2)  to  (A6) 
for  the  (constant)  temperature,  humidity,  and  pressure  appropriate  for  each 
path  segment.  The  path  length  segments  are  assumed  to  be  short  enough  that 
the  calculated  absorption  coefficient  Is  a good  representation  of  the  average 
absorption  over  the  segment. 

For  sound  propagation  from  the  same  source  and  over  the  same  path,  the 

squared  sound  pressure  at  the  receiver  after  propagation  through  atmosphere 
2 

2,  Pj^2’  written  as 

p‘2  (f.C)  = [pg(f,0)](e  •^)(e  ^)  •••  (e  ") . (A14) 

Substituting  equations  (A13)  and  (A14)  In  equation  (A12)  gives 
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TR 


2 (“2,r“2.2^  ^^2 

2 ^“n.r“n.2> 


(A15) 


"4.TR  = 10  1°8  ^ 


(A16) 


where  A is  the  absorption  factor  defined  by  the  exponential  term  on  the  right 
hand  side  of  equation  (A15). 

To  apply  equation  (A15),  the  absorption  coefficients  a and  a are 

O ) X XX  f ^ 

calculated  using  subroutine  ANSAB  for  each  pathlength  segment  AC  for 

n 

specified  values  of  temperature,  humidity  and  pressure  appropriate  to  that 
segment.  The  separate  terms  in  the  exponent  are  calculated  and  summed  to 
determine  the  value  of  A and  hence  AL^  XR*' 


An  alternate  expression  for  the  absorption  factor  in  equation  (A15) , 
which  is  more  useful  for  some  applications,  can  be  obtained  from  the  mathe- 


matical identities 


-2a  AC  -(a  /lO)  AC 

n ^ /is  n n 

e =10 


(A17) 


2a  AC  (a  /lO)  AC 

n n n n 

e =10 


(A18) 


With  equations  (A17)  and  (A18),  the  absorption  factor  becomes 


A = 10 


+ (a„  ,-a_  ,) to/10)) 
n,  i n,  Z n 


(A19) 


*In  applying  these  pure-tf ne  atmospheric  absorption  adjustment  factors  to 
1/ 3-octave-band  sound  pressure  levels  containing  discrete-frequency  com- 
ponents, it  is  assumed  that  the  tone  frequency  is  always  within  the 
essentially  zero-loss  portion  of  a filter's  passband. 
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ATMOSPHERIC  ABSORPTION  ADJUSTMENT  FOR  BANDS  OF  NOISE 


Analysis 

Reference  1 contains  a discussion  of  the  application  of  the  proposed  new 
technique  of  ANS  SI. 26  for  determining  pure-tone  sound  absorption  coefficients 
to  calculation  of  the  atmospheric  absorption  losses  for  noise  analyzed  by 
constant-percentage-bandwidth  or  fractions 1-octave-band  filters.  The  procedure 
in  reference  1 for  determining  band-loss  adjustment  factors  is  not  directly 
applicable  to  the  problem  of  correcting  for  differences  in  sound  absorption 
losses  occurring  during  propagation  through  a horizontally  stratified  atmos- 
phere. The  analysis  described  in  this  Appendix  uses  some  of  the  methods  pre- 
sented in  reference  1 to  develop  a procedure  applicable  to  adjustment  of 
measured  aircraft  noise  levels  from  test  to  reference  meteorological  conditions. 
The  basic  concepts  for  the  analysis  procedure  for  bands  of  noise  were  developed 
by  Louis  C.  Sutherland  for  the  Sl-57  Working  Group  that  prepared  the  proposed 
new  American  National  Standard  on  atmospheric  absorption. 


The  sound  source  Is  assumed  to  produce  a noise  signal  having  a continuous 
pressure  spectral  density  over  the  complete  range  of  frequencies  Included  in 
the  analysis.  The  pressure  spectral  density  of  the  noise  signal  at  the  source 
is  denoted  by  Gg(f,0)  at  any  frequency  f.  After  propagation  from  the  source 
to  the  receiver  over  distance  C,  the  pressure  spectral  density  of  the  noise 
signal  is  G^(f,C)  for  propagation  through  atmosphere  1 with  test  meteorological 
conditions  and  for  propagation  through  atmosphere  2 with  reference 

meteorological  conditions. 


By  analogy  to  the  development  of  equations  (A13)  and  (A14)  in  Appendix  A, 
the  pressure  spectral  densities  G^^Cf,?;)  and  G2(f,?)  can  be  related  to  Gg(f,0) 
by 


Gg(f,0)  (e 


■2“l,l  ^^1 


) (e 


"^“2,1  ^^2. 


-2a  . AC 

(e  ") 


(Bl) 


and 


G2(f,C)  = Gg(f,0)  (e 


-2«1,2  ^^1 


) (e 


-2“2,2  ^^2, 


-2a  - AC 

(e  (B2) 
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We  relate  G to  by  solving  equation  (Bl)  for  G , substituting  the 
2 1 ^ 

result  into  equation  (B2) , and  combining  related  terms  in  the  exponent  to 
obtain 

[2(a  - a „)  AC  + 2(a  - a„  ) AC,  + ••• 

G^Cf.C)  = G^(f.C)  {e 


+ 2 (a 


n,l 


-a  ) AC  ], 
n, 2 n > 


(B3) 


= G^(f,C)  A(f,C)  (B4) 

where  we  have  used  the  same  symbol  A for  the  absorption  factor  defined  pre- 
viously by  equations  (A15)  and  (A16)  for  pure- tone  sound  signals. 


The  set  of  1/3-octave-band  sound  pressure  levels  after  propagation  through 
atmosphere  1 and  analysis  by  filters  having  response  characteristics  1 is 
called 

The  set  of  1/ 3-octave-band  sound  pressure  levels  at  the  receiver  location 
after  propagation  through  atmosphere  2 and  analysis  by  1/ 3-octave-band  filters 
having  transmission  response  2 is  called  L2-  We  want  to  find  L2  knowing  the 
values  of  L^,  the  characteristics  of  the  atmosphere,  and  the  characteristics 
of  the  filters. 


The  measurement  process  that  results  in  the  measured  levels  and  that 
would  result  in  the  levels  L2  is  shown  schematically  in  figure  Bl.  The  sound 
pressure  signal  incident  on  the  microphone  is  a broadband  signal  whose  amplitude 
p^(t)  varies  randomly  with  time.  The  frequency  spectrum  of  the  signal  Pj^(t) 
is  a smooth,  continuous  function  as  indicated  by  the  sketch  of  Gj^(f)  ^ f. 


The  sound  pressure  signal  is  recorded  and  processed  by  some  data-acqulsl- 

tion/data-processing  system.  The  system  contains  a set  of  1/3-octave-band 

2 

filters  each  of  which  has  a power  transfer  function  |H^(f;f^)|  defined  by  the 
square  of  the  absolute  value  of  its  frequency  response.  The  response  curve  for 
one  of  the  filters  with  band  center  frequency  f^  is  shown.  Except  for  the 
filter,  the  rest  of  the  entire  data-acquisitlon/data-processing  system  is 
assumed  to  have  unity  gain  at  all  frequencies  of  interest. 
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Figure  B1.- Illustration  of  steps  involved  in  determining  1 /3-octave-band  sound  pressure  levels  representing 
the  spectrum  of  a broadband  random  noise  incident  on  a microphone. 
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Depending  on  the  bandwidth  of  the  filter,  the  signal  at  the  output  of  the 
filter  has  an  approximately  sinusoidal  waveshape  but  with  a random  amplitude  as 
indicated  by  the  sketch  of  p(t)  t.  The  spectral  density  of  the  voltage 
equivalent  of  the  p(t)  signal  is  shown  by  the  sketch  of  S]^(f;f^)  va  f around 
the  band  center  frequency  f^.  The  function  Sj^(f;f^)  is  related  to  the  pressure 
spectral  density  of  the  input  signal  G^(f),  after  conversion  from  a pressure 
signal  to  a voltage  signal  and  passage  through  a filter,  by 

S^(f;f^)  = G^(f)  |H^(f;f^)|2.  (B5) 

For  the  example  shown  in  figure  Bl.the  Sj^(f;f^)  function  is  equal  to 
G^(f)  in  the  flat-response  (passband)  region  of  the  filter  around  the  center 
frequency  f^.  In  the  frequency  region  below  the  lower  bandedge  frequency  f^^ 
(i.e.,  in  the  lower  stopband),  the  S^(f;f^)  function  decreases  more  slowly 
than  it  does  in  the  frequency  region  above  the  upper  bandedge  frequency  fy 
(i.e.,  in  the  upper  stopband)  because  the  Gj^(f)  function  increases  at 
frequencies  less  than  fy  and  decreases  at  frequencies  greater  than  fy.  Unless 
the  rolloff  rate  of  |Hi(f;f^)|  is  steeper  than  the  slope  of  Gj^(f),  the 
Si(f;fi)  function  will  be  assymmetric  and,  depending  on  the  shape  and  slope 
of  G^(f),  there  can  be  significant  contributions  from  either  the  lower  or 
upper  stopband  frequency  regions.  If  the  filter  were  an  ideal  filter  with 
unity  gain  in  the  passband  and  infinite  rejection  in  the  stopbands,  then 
Si(f;fi)  would  just  equal  G^(f)  between  the  passband  limiting  frequencies  f^ 
and  fy  and  be  zero  at  all  other  frequencies. 

The  final  step  in  the  process  is  the  detection  of  the  equivalent  p(t) 

signal  and  the  display  (or  output)  of  the  sound  pressure  level  for  the  band. 

Detection  involves  squaring  and  time  averaging  to  find  the  mean-square  (or 

2 

root-mean-square)  value  of  p(t)  or  < p (t)  > where  the  <>  signs  indicate  time 
averaging  and  it  is  understood  that  suitable  conversion  or  scaling  factors 
would  be  included  to  provide  proper  units. 

The  sound  pressure  level  in  band  at  f^  for  the  filtered 

equivalent  sound-pressure  signal  p(t)  is 
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» 10  log  [•^P^>/PrefJ 

2 

where  P^g£  Is  the  square  of  the  reference  pressure  “ 20  pPa). 

The  mean-square  value  of  a signal  is  equal  to  the  integral  over  frequency 
of  the  pressure  spectral  density  of  the  signal.  Thus,  equation  (B6)  can  be 
written  as 

L^(f^)  = 10  log  {[Q/"Sj(f;f^)df]/pJg^}  (B7) 

or,  with  equation  (B5),  as 

Lj(f^)  - 10  log  {[Q/"[G^(f)][|Hj^(f;f^)|^]df]/p2^^}.  (B8) 

Similarly,  for  the  sound  pressure  levels  L2(f^)  after  propagation  through 
the  reference  atmosphere,  we  can  write 

, ^ L2(fj)  - 10  log  {[Q/”s2(f;f^)df]/p^^^}  fB9) 

or 

L2(f^)  • 10  log  {[Q/“[G2(f)](|H2(f;f^)|^]df]/p^^^}.  (BIO) 

As  In  the  Initial  statement  of  the  problem,  we  assume  that  the  same  set  of 
filters  would  be  used  to  measure  both  and  L2  and  thus  set 

|Hj^(f;f^)|^  - |H2(f;f^)t^-  |H(f;f^)|^.  (Bll) 

We  also  use  equation  (B4)  to  replace  G2(f)  In  equation  (BIO)  to 

obtain 

I 

L2(f^)  - 10  log  {[Q/“G^(f)  A(f)  |H(f;f^)|^df]/p^^^}  (B12) 

and 

L^(fj)  - 10  log  {[Q/“G^(f)  |H(f;f^)l^df]/p^^^}.  (B13) 

The  band-loss  atmospheric-absorption  adjustment  factor  from  test  to 

reference  meteorological  conditions,  AL„  „,  is  the  quantity  we  seek.  It  is 

Oy  XK 

determined  for  the  band  with  center  frequency  f^^  from 

%TR  ■ "2  <'l>  - h ^'l> 


J 
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and  Is  to  be  added  to  the  measured  levels  L^(f^)  to  obtain  the  adjusted 
levels  L2(f  . 

For  any  band,  using  equations  (B12)  and  (B13),  we  can  write 

ALg  ‘ 10  log  {[q/^  a |h|^  df]/[Q/”Gj  |h|^  df]}  (B15) 

where  the  notation  showing  the  functional  dependence  on  frequency  has  been 
omitted. 

We  now  assume  that  the  slope  of  the  pressure  spectral  density  is  not 
'too  steep'  (neither  too  positive  nor  too  negative)  and  that  the  response 
|h|  of  each  of  the  real  filters  can  be  adequately  represented  by  that  of  ideal 
filters  with  the  same  set  of  band  center  frequencies.  In  other  words,  we 
assume  that  there  is  negligible  contribution  to  the  integrals  in  equation  (B15) 
at  frequencies  outside  the  filter  passband  and  replace  the  infinite  frequency 
range  by  an  integration  from  f to  f with  |h|  = 1.0. 

Determination  of  the  range  of  pressure  spectral  density  slopes  for  which 
the  response  curves  of  practical  filters  could  be  safely  approximated  by  the 
response  of  an  ideal  filter  was  beyond  the  scope  of  this  project.  This  matter 
should  be  investigated  in  actual  applications  of  the  techniques  discussed  here 
because  sound  spectra  with  large  positive  or  large  negative  slopes  can  result 
in  significant  energy  being  transmitted  by  the  stopbands  of  certain  practical 
filters.  Large  negative  slopes  are  often  encountered,  for  example,  in  the 
high-frequency  part  of  aircraft  noise  spectra  at  long  distances  because  of 
atmospheric  absorption  effects. 

With  the  assumptions  described  above,  equation  (B15)  becomes 

ALg,TR  = 10  log  {[^f  Gj^  A df]/[^/  . df]}.  (B16) 

L L 

To  Integrate  equation  (B16)  we  require  an  expression  fcr  the  pressure 
spectral  density  function  Gj^(f).  We  assume  that  we  can  find  a reasonable 
approximation  to  the  true  pressure  spectral  density  function,  over  the 
frequency  range  of  a filter's  passband,  by  means  of  a line  segment  of  slope 
£'  as  shown  in  the  sketch  in  figure  B2.  At  the  band  center  frequency  f^^,  the 
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Figure  B2.-Exatnple  of  approximation  to  , ^assure  tpectral  deiMity  G^(fi  by  litte 

segment  of  slope  i ' over  frequency  range  from  f to  fg  of  a 1/3-octave- 
band  filter  with  band  center  frequency  fj. 

line  segment  approximation  has  the  value 

Along  the  line  segment  in  figure  B2,  the  slope  between  Gj^(f)  at  any 
general  frequency  f and  at  band  center  frequency  f^  can  be  written  as 

A'  - [log  Gj^(f)  - log  G^(f^)]/[log  f - log  f^].  (B17) 

From  equation  (B17)  we  can  write  the  general  expression  for  the  straight- 
line  approximation  to  Gj^(f)  as 


G^(f)  - G^(fpt(f/f^)^'l 

(B18) 

G^(f)  - G^(f^)  f"*’  f^'. 

(B19) 

Substituting  equation  (B19)  into  equation  (B16)  yields,  after  pulling  the 
-I' 

constant  ^^^^l^^l  outside  the  integrals  and  canceling. 


- 10  log  {[^f  ^Af^’  df]/[^A  f^’  df]} 
L L 


(B20) 
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I 


t 


To  evaluate  equation  (B20)  we  begin  with  the  denominator  term.  When 
I'  ft  - 1,  we  have 


/ ^ f^'df  = [l/()l'+l)] 


When  i'  = -1,  we  have 


/ ^ f"^df  = In 
f 

L 


where  In  represents  natural  logarithms. 


(B21) 


(B22) 


To  integrate  the  numerator  term  in  equation  (B20)  we  replace  the  single 
Integral  over  the  frequency  range  of  the  passband  by  a summation  of  integrals 
over  a number  of  equally-spaced  constant-percentage-bandwidth  segments  SE. 

This  method  of  integration,  instead  of  a straightforward  numerical  integration 
technique  available  from  a number  of  standard  scientific  computational 
subroutines,  was  selected  to  minimize  computer  time  (and  cost)  for  a given 
level  of  accuracy. 


Thus,  we  write  for  the  numerator  term 

/ “ A (’•'  df  - .E  ,/  A f*'  d£.  (B23) 

The  number  of  segments  SE  is  chosen  to  be  large  enough  that  the  frequency 
spectrum  of  the  absorption,  over  the  length  A?  of  each  step  along  the  propaga- 
tion path,  can  be  assumed  to  be  represented  by  a series  of  elements  over  which 
the  absorption  is  proportional  to  frequency  raised  to  some  power  k.  By  absorp- 
tion over  a pathlength  AC,  we  mean  the  ratio  of  the  squared  sound  pressures 
at  the  beginning  and  end  of  the  step,  l.e.,  equations  (A17)  or  (A18) . The 
absorption  spectrum  is  most  conveniently  represented  using  equation  (A18)  on 
logarithmic  scales  and,  over  each  of  the  narrow  bands  of  frequency  in  the 
filter  passband  defined  by  the  number  of  segments  SE,  the  absorption  spectrum 
is  approximated  by  a straight  line  with  slope  k. 
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Thus,  for  a given  temperature,  humidity,  pressure,  and  pathlength  distance, 
the  slope  k of  the  straight-line  approximation  from  a base  or  reference  fre- 
quency fj  to  any  general  frequ'“.ncy  f > fj  can  be  written  as 

a(f)AC/lO 


a(f  )AC/10 

k = [log  10  - log  10  ^ ]/[log  (f)  - log  (fj)]  (B24) 


or 


k = {[a(f)  - a(fj)](AC/10}/log  (f/fj) 


(B25) 


A general  expression  for  any  term  in  the  absorption  factor  A defined  by 
equation  (A19)  can  then  be  written,  from  equation  (B24),  as 


a(f)AC/10  a(f  )AC/10 

10  = 10  ^ 


(f/fj) 


(B26) 


Substituting  the  general  expression  from  equation  (B26)  into  equation 
(A19)  for  each  of  the  separate  meteorological  conditions  and  pathlength 
segments  yields 


A - 10 


+ •••  + [a  i(f.)  - a ,(f,)]  (AE  /lO)} 
n,l  j n,2  j n 


[(k  . - k ,)  + (k,  , - k ) + •••  + (k_  - k_  „)] 


1,1  1,2  2,1  2,2 


n,l  n,2' 


(B27) 


We  Introduce  the  abbreviations 


A(fj)  - 10 


+ •••  + [a„  ,(f,)  - a^  ,(f,)](A£  /lO)}  (B28) 

n,i  J n,2  j n 

for  the  value  of  the  effective  absorption  factor  at  aad 

K - - k^^2)  + (^^2,1  " *^2,2^  + •••  + (k^^^  - k^^^)  (829) 

for  the  effective  slope  parameter.  The  terms  A(fj)  and  K represent  effective 
or  average  values  over  the  total  propagation  path  resulting  from  differences  in 
atmospheric  absorption  between  test  and  reference  conditions. 
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With  equations  (B28)  and  (B29) , equation  (B27)  becomes 


A - A(f^)  [(f/f^)n, 


(B30) 


We  can  test  the  reasonableness  of  equation  (B30)  by  noting  that  If  the 
test  and  reference  meteorological  conditions  were  Identical,  the  value  of 
A(fj)  would  be  1.0  and  K would  be  0.0  giving  A = 1.0  as  It  should. 

With  equation  (B30) , one  of  the  Integrals  In  the  sum  on  the  right  hand 
side  of  equation  (B23)  becomes 


or 


A(fj)  [(f/fj)*^]  f^’  df 

A(t,)  £■'«'  df. 


(B31) 


(B32) 


When  K+H'  ^ -1,  we  have 

A(£.)  ff  £'‘«’d£  . A(f.)  £-''  [l/(K+l'+l)l  (f™'”  - ff'’*!) 

3 3 3 j 3"^*  J 

= A(fj)  [l/(K+£'+l)]  -1] 

= A(fj)  [l/(K+Jl’+l)][(fj^^/fj)*^'^^''^^  -1].  (B33) 

When  K+i.'  = -1,  we  have 


A(fj)  f"*^  f"^  df  - A(f^)  f'*^  In  (fj^j/f^). 


(B34) 


With  equations  (B21) , (B22) , (B33),  and  (B34)  there  are  thus  four  possible 
cases  for  the  ratio  of  Integrals  In  equation  (B20). 


58 


APPENDIX  B 


When  I’  ^ -1  and  K+£'  ^ -1,  equations  (B33)  and  (B21)  yield  the  ratio 

£.'+1  K+t  '+1 

Z A(f  ) r - 1] 

i=l  J J J . ( 

[l/(£'+l)] 


(B35) 


For  convenience,  we  divide  numerator  and  denominator  of  equation  (B35) 
£ '+1 

by  and  move  the  constant  terms  from  the  denominator  to  under  the 

summation.  This  manipulation  gives 
SE 

E A(fj)[(fj/f^)*'''^^]  [(£'+!)/ (K+£'  + l)][(fy/f^)^’'^^-(fj^/f^)^’‘^^]"^ 


V’J.tf  » +1 

- 1]  = bli^ 


"S.TR  = 10 


(B36) 


(B37) 


when  £'  ^ -1  and  K+£'  ^ -1,  or  K - (£'+l). 

For  the  second  case,  we  take  £’  ^ -1  and  K+£'  = -1,  i.e.,  K = - (£’+l). 
Equations  (B34)  and  (B21)  yield 


E A<fj)  tj-  1„ 

JlZ± I 

o*+i  P,*+1 

[l/(£’+l)]  (fj  - f‘ 

Employing  the  same  technique  as  before  and  using  K = - (£'+l)  gives 


(B38) 


£’+1,-1 


E A(fj)  [(fj/f^)"  (£'+l)  [(Vfp  - (V^i^  J 1" 


SE 

E BL2 


(B39) 
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"S.TR  = 

3=1 

when  V ^ -1  and  K+Jl'  = -1,  or  K = - (£'+l). 


(B40) 


For  the  third  case,  we  have  V = -1  and  K+i’  -1  or  K 0.  Equations 
(B33)  and  (B22)  yield 


I A(f  ) (1/K)  J - 1] 

3 = 1 ^ ^ ^ 


<B41) 


KJ  V 

E A(f.){l/[K  ln(f,,/f,)]}[(f,.,/f,)'^  - 1]  = E BL3 
j=l  J ^ ^ J J j = l J 


(B42) 


"S.TR  = BL3j 

when  I'  = -1  and  K+H'  ^ -1  (or  K 0). 


(B43) 


The  final  case  is  when  = -1  and  K+Jl'  = -1,  i.e.,  when  K = 0.  With 


K = 0,  we  have 


"4.TR  = O-O- 


(B44) 


Note  that  when  K=0  the  adjustment  factor  is  0.0  regardless  of  the  noise 

SE 

slope  For  -1  and  K=0  equation  (B36)  yields  E BLli  = 1.0  and  hence 

j=l 

^Lb.TR  “ 0.0  by  equation  (B37);  for  £=  -1  and  K=0,  equation  (B44)  applies. 

To  simplify  equations  (B36),  (B39),  and  (B42)  we  note  that  the  ratio  of 
upper  to  lower  bandedge  frequencies,  is,  by  definition,  a constant  for 

any  constant-percentage-bandwidth  filter.  We  denote  this  frequency  ratio  by 
the  symbol  RF. 


The  upper  bandedge  frequency  is  related  to  the  band  center  frequency  by 


fu  = f,  RF^ 


(B45) 
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The  lower  bandedge  frequency  is  related  to  the  band  center  frequency  by 

Since  the  passband  is  divided  into  SE  equally  spaced  frequencies  (on  a log 
frequency  scale),  the  ratio  of  any  two  consecutive  frequencies  is 

fj_j^^/fj  = (B47) 

The  choice  of  the  value  of  the  SE  should,  ideally,  consider  the  slopes 
of  the  noise  spectrum  and  the  absorption  spectrum.  However,  on  the  basis  of 
experience,  subdivision  of  the  1/3-octave  bands  into  twenty  constant-percentage 
intervals  was  estimated  to  be  able  to  provide  adequate  accuracy  for  practical 
cases  of  Interest  to  aircraft  flyover  noise  measurements. 


In  order  to  provide  a set  of  regularly  spaced  band  center  frequencies 
that  do  not  allow  irregular  overlapping  of  the  frequency  response  of  the 
contiguous  filters,  the  band  center  frequencies  were  defined  to  be  the 
precise  frequencies  for  1/3-octave-band  filters,  as  specified  in  reference  12, 
namely 

f.  . 10™'”  - 10"'“  (B48) 

where  N is  any  integer:  positive,  negative,  or  zero.  The  range  of  N for 
these  1/3-octave-band  analyses  is  from  17  to  40  covering  nominal  band  center 
frequencies  from  50  to  10,000  Hz. 

With  the  precise  band  center  frequency  defined  by  equation  (B48)  the 
band  frequency  ratio  is 

RF  = = lO^'^^®  = 1.2589.  (B49) 

Note  that  equations  (B48)  and  (B49)  define  a set  of  frequencies  that  coincide 
exactly  with  the  internationally  standardized  and  preferred  band  center 
frequencies  at  the  Important  decade  frequencies  of  100,  1000,  and  10,000  Hz. 


Once  a value  has  been  established  for  RF,  then  the  frequencies  fj  depend 
only  on  the  choice  of  the  number  of  segments  SE,  and,  as  stated  earlier. 
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acceptable  accuracy  can  be  obtained  with 


SE  = 20. 


1/SE  1/20 

The  size  of  the  frequency  steps  between  f . and  f ^ Is  RF  = 10 

1/3  ^ 

1,0116  for  SE  = 20.  [If  RF  were  2 , then  the  frequency  Interval  would  be 

2I/6O  _ 1/60-th  octave  bands  for  SE  = 20.] 


(B50) 

1/200 


The  set  of  fj  frequencies  begin  at  f^^  and  Increment  by  the  frequency 
spacing  factor  RF^^^^,  Thus,  the  f^  frequencies  are  specified  by 

f /f^  = (RF"^'^^) 


(B51) 


When  j = 1,  f./f,  = iJf..  When  j = SE,  f./f.  = f„p/f,  = (RF 
jl  ^^1/2  SEl 

and  f... /f.  = fop.i/fj  = RF  = f,,/f.  at  the  end  of  the  Interval  for  the 
j+1  1 bti+i  1 U 1 

passband. 

The  only  remaining  Item  needed  to  determine  the  band-loss  adjustment  factors 
Is  specification  of  the  spectral  slope  V . Recall  that  £■'  is  the  slope  of  the 
pressure  spectral  density  of  the  sound  pressure  at  the  microphone  before  being 
analyzed  by  the  1/3-octave-band  filters.  This  slope  Is  not  available  since 
the  pressure  spectrum  of  the  source  sound  pressure  Is  not  known;  only  the 
filtered  band  levels  are  available. 

To  proceed,  we  assume  that  the  slope  of  a straight  line  between  the  average 
pressure  spectral  densities  W^^  in  the  bands  below  and  abo/e  the  band  of 
interest  at  center  frequency  f^  is  a reasonable  approximation  to  the  slope  of 
the  true  pressure  spectral  density  G^.  The  average  pressure  spectral  density 
over  the  passband  of  the  filter  at  center  frequency  f^  is 


W^(fl)  = <p^(f^)>/BW^ 


(B52) 


where  <p^(f^)>  is  the  mean-squared  sound  pressure  in  the  band  at  f^  and  BW^ 

is  the  bandwidth  of  the  band  or  (f  -f  ).  With  equations  (B45)  and  (B46)  and 

U L> 

equation  (B6),  we  can  express  equation  (B52)  as 


2 L (f,)/10  -1/2 

W^(fi)  = 10  ]/[f^(RF^'^  - RF 


(B53) 
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Similarly, 


and 


j 1/2  -1/2 

Wi^^i+1^  “ ^Pref  l/tfi+i(RF  ' - RF  (B54) 

2 1/2  -1/2 

Wi(fi-i)  = 10  ^ ^ ^ ]/[f^_;^(RF  ^ - RF  ' )].  (B55) 


The  slope  i'  of  the  straight  line  over  the  band  at  between 
at  '^l^^i  1^  ^i-1  shewn  by  the  sketch  in  figure  B3.  This  slope 

is  determined  from 


i'  = [log  - log  W^(f^_j^)]/[log  f^^^  - log  f^_^].  (B56) 

With  equations  (B54)  and  (B55),  equation  (B56)  becomes 

I'  = [{[L^(f^^P  - L^(f^_^)]/10}  - log  (f^^^/f^_^)]/log 
a'  = [{[L^(f^^^)  - L^(f^_j^)]/10}/log  - 1-  (B57) 


LOG  f 

Figure  B3.-Example  of  approximate  determination  of  slope  /'from  average  pressure 
spectral  densities  in  bands  below  and  above  the'  band  of  interest  at 
center  frequency  f j. 
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The  slope  of  measured  sound  pressure  over  the  band  at  from  the  band 
above  to  the  band  below,  is  the  band  noise  slope  on  a plot  of  the  logarithm 
of  the  squared  sound  pressure  versus  the  logarithm  of  band  center  frequency. 

Thus  the  band  noise  slope  over  the  band  at  f^  Is 

A = {[L^(f^^^)  - L^(f^_^)]/10}/log  (B58) 

where  the  factor  of  10  Is  Included  to  eliminate  the  decibels  and  provide  a 
measure  of  the  mean-square  sound  pressure  In  the  band. 

With  equation  (B58) , the  spectral  slope  I'  can  be  written  as 

1'  - 1 - 1.  (B59) 

With  RF  defining  the  frequency  ratio  for  any  set  of  band-spaced  frequencies, 

2 

we  have  RF  and  = f^/RF.  Thus  ~ ^ write 

the  band  noise  slope  of  equation  (B58)  as 

a = - L^(f^_j^)]/20  log  (RF)  (B60) 

for  every  band  except  for  first  and  last  band. 

The  slope  Z over  the  first  band  Is  found  from  the  difference  In  sound 
pressure  level  between  that  of  the  second  band  and  that  of  the  first  band. 

The  slope  over  the  last  band  Is  found  from  the  difference  In  sound  pressure 
level  between  that  of  the  last  band  and  that  of  the  next-to-last  band. 

Thus  for  the  first  band 

Z = - 4(^1)  1/10  log  (RF)  (B61) 

and  for  the  last  band 

Z = [^(f^)  - Lj^(fi_i)]/10  log  (RF).  (B62) 

Equations  (B59),  (B60),  (B61),  and  (B62)  define  the  spectral  slopes 
Z and  Z'  . 

With  equations  (B59),  (B45) , (B46),  and  (B47) , the  band  loss  adjustment 
factors  can  be  written  In  a more  convenient  form  as  follows. 
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When  jt  -1,  [£  ^ 0],  and  K + -1,  [(K+Jl)  ^ 0 or  K ^ -I,  but 

K 0],  we  have  from  equation  (B36) 

BLlj  - A(fj)[(fj/f^)^][i/(K+«,)][l/(RF*''^  - RF"^'^^)]  - 1].  (863) 

When  1'  ^ -1,  [<,  ^ 0],  and  K + il'  - -1,  [K  + )l  - 0 or  K - -t,  but  K ^ 0], 
we  have  from  equation  (839) 

BL2j  = A(fj)[(fj/f^)*']  (i/SE)  [1/(RF*'^^  - RF"*'^^)  ] In  (RF).  (B64) 

When  I'  = -1,  [il=0],  and  K+J!.'?^  -1,  [(K+il?*0  or  Kj*0],  we  have  from 
equation  (B42) 

BL3j  = A(fj)  {1/[K  In  (RF)  ] } (RF*^^®®^  - 1).  (B65) 

Since,  by  the  definitions  we  have  selected,  the  spectral  slopes  I'  or  1 
are  constant  over  the  passband  of  the  filter,  the  three  cases  described  by 
equations  (B63),  (B64),  and  (B65)  really  only  define  two  options,  assuming 
K?^0  at  any  frequency  f^  in  any  passband.  If  K“0,  then  there  is  no  adjustment 
factor  and  equation  (B44)  applies  regardless  of  noise  slope. 

The  two  cases  are  for  when  l<^0  and  when  fc“0.  When  1^0,  then  either 
equation  (B63)  or  equation  (B64)  holds  depending  on  whether  K+J!,?to  or  K+)l“0  at 
a particular  frequency  step  In  the  summation  over  a passband.  We  therefore 
define  an  Index  counter  BLJl  to  store  contributions  to  the  summation  when 
if^O  and  where 


BLJl  = BLlj  when  K+Hi^O,  eq.  (B63), 
and  BLJl  « BL2j  when  K+1-0,  eq.  (B64). 

The  band  loss  adjustment  factor  when  is  then 


SE 

ALg,TR  - 10  log  E BLJl. 


(B66) 

(B67) 


(B68) 
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When  £ = 0,  and  K+il?t0  or  Kj^O,  then  only  equation  (B65)  applies  and  no 
additional  storage  counter  Is  needed.  Thus 

SE 

aL_  = 10  log  E BL3,  (B69) 

B,TK  3 

when  1=0  and  K?*0. 


Recapitulation 

The  various  expressions  needed  to  carry  out  the  calculations  of  band 
loss  adjustment  factors  are  repeated  here  for  convenience. 

When  the  band  noise  slope  1 is  nonzero  (1#0),  then  the  band-loss  adjust- 
ment factor  AL„  „ is  found  from 
Dy  IK 


SE 

AL_  _ = 10  log  E BLJl 
B.TR 


where 


when  K+lj^O,  or  from 


BLJl  = BLl 


i 


BLJl  = BL2 


j 


(B68) 


(B66) 


(B67) 


when  K+1  = 0,  and  where  K is  the  effective  sound-absorption-factor  slope  over 
the  entire  sound  propagation  path. 


When  1=0  and  K?*0,  then 


SE 


When  K=0 , then 


"S.TR  = «"3j' 


"S.TR  “ 0-0 


(B69) 


regardless  of  the  band  noise  slope  1. 
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The  terms  in  the  summations  are  defined  by 


BLlj  - A(fj)[(fj/f^)^]()l/(K+l)][l/(RF^/2  - RF  <^/2)  j - 1]  (B63) 


and 


BL2j  - A(fj)[(fj/f^)^](£/SE)[l/(RF^^^  - RF"*"^^)]  In  (RF)  (B64) 


and 


BL3j  = A(fj){l/[K  In  (RF) ] } (RF 


,K/SE 


1).  (B65) 

The  various  expressions  in  the  three  sum  terms  are  defined  as  follows. 


The  absorption  factor  A(f^)  at  one  of  the  frequencies  f^  in  the  numerical 


integration  over  the  frequency  range  of  the  filter  passband  from  the  lower  f 


to  the  upper  fy  bandedge  frequency  is 


A(fj)  = 10 


{[ai,l(fj)  “ (AC^/IO)  + (AC2/IO) 


+ •••  + [a^  -(f.)  - a„  ,(f,)]  (At  /lO)} 
n,i  j n,2  j n 


(B28) 


where  the  various  values  of  absorption  coefficient  a are  determined  by 


equations  (A2)  to  (A6)  at  frequency  fj  for  the  temperature,  humidit  , and 


pressure  conditions  appropriate  for  each  segment  along  the  length  of  the 
propagation  path  (segments  of  length  AC)  for  either  the  test  atmosphere 
(second  subscript  1)  or  reference  atmosphere  (second  subscript  2). 


The  effective  absorption-factor  slope  term  K,  over  the  frequency  range 
from  fj  to  aiid  over  the  entire  length  of  the  propagation  path,  is 


K = 2^  ^'^2  1 ~ ^^2  2^  ^ 


+ (k  - k ) 
n,l  n,z 


(B29) 


while  the  individual  absorption  slopes  k are  determined  at  the  appropriate 
temperature,  humidity,  and  pressure  from 


a/SE, 


k - {[a(fj^^)  - a (fj)]  (AC/10}/log  (RF'"''""') 


(B70) 


after  subsituting  for  f in  equation  (B25)  and  making  use  of  equation 


(B47)  for  the  frequency  ratio  f^^^^/fj. 
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The  1/ 3-octave-band  frequency  ratio  is 


RF  = 10 


(B49) 


and  the  number  of  segments  over  the  frequency  range  of  the  filter  passband 


SE  = 20. 


(B50) 


The  Initial  value  of  the  fj  frequencies  is  the  lower  bandedge  frequency 


fj^  = = 10 


(B46) 


where  f^  Is  the  band  frequency  defined  by' 

f^  = 10*^^^®  (B48) 

and  where  N ranges  from  17  to  40  for  the  24  bands  considered  here. 

The  final  value  of  the  f^  frequencies  is  the  upper  bandedge  frequency 


fy  = RF^'^^f^  = lO^'^^^f^. 


(B45) 


The  Increment  in  f^  frequencies  is  determined  from 


(B47) 


while  a general  f^  frequency  is 


fj  = f^  (RF"^*^^) 


(B51) 


For  every  band  except  the  first  and  last  bands,  the  band  noise  slope  H 


is  found  from 


a = - Lj<fi-i)l/20  log  (RF), 


For  the  first  band,  we  use 


Jl  = [L^(f2)  - L(f^)l/10  log  (RF) 


(B60) 


(B61) 
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and  for  the  last  band  we  use 

il  = [Lj(f^)  - l°g 

The  calculation  procedures,  symbols,  and  terminology  given  here  In  this 
recapitulation  plus  equations  (A2)  to  (A6)  from  Appendix  A are  used  directly 
for  the  computer  program  BANDLOSS,  see  Appendices  C and  D. 
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LISTING  OF  COMPUTER  PROGRAM  SOURCE  STATEMENTS 


Source  statements,  in  FORTRAN  IV,  for  program  BANDLOSS  are  listed  in  this 

Appendix.  The  program  was  written  to  calculate  atmospheric-absorption 

adjustment  factors  for  a single  0.5-second  sample  of  an  actual  aircraft 

flyover  noise  measurement.  The  aircraft  noise  spectrum  is  assumed  to  be  that 

of  a continuous  broadband  free-field  signal  with  no  discrete-frequency  pure- 

tone  components.  The  program  determines  1/3-octave-band  absorption  adjustment 

factors  AL„  Subroutine  ANSAB  is  called  to  calculate  atmospheric 

Oy  IK 

absorption  at  specific  frequencies  within  the  passband  of  each  1/3-octave 
filter  for  each  segment  of  the  propagation  path. 

Input  data  are  either  defined  in  a data  list  or  are  read  in.  The  data 
list  defines  the  nominal  center  frequencies  of  the  1/3-octave  bands  from  50  Hz 
to  10,000  Hz  for  use  in  labelling  the  output  data. 

Input  data  are  read  in  either  directly  from  the  terminal  keyboard  or  from 
input  files  (if  operating  online).  They  could  also  be  read  in  from  punched 
cards  or  they  could  be  supplied  from  another  computer  program  which  incorporat- 
ed BANDLOSS  as  a subprogram.  The  terminal  input  data  consist  of:  IMAX,  TS, 

TN,  TOH,  AH,  ZM,  AS,  and  the  set  of  24  1/ 3-octave-band  sound  pressure  levels, 
SPL.  Measured  meteorological  data  applicable  to  the  input  SPLs  are  read  from 
the  defined  files  as  TC,  HR,  and  BP. 

Output  data  consist  of  tabulated  values  of  the  input  1/ 3-octave-band 

sound  pressure  levels,  SPL,  the  adjusted  1/ 3-octave-band  sound  pressure 

levels,  SPLA,  and  the  1/ 3-octave-band  atmospheric-absorption  adjustment 

factors,  BLTR  = AL  . Associated  airplane  and  geometric  parameters  are 
B f TK 

listed  with  SPL,  SPLA,  and  BLTR.  Meteorological  parameters  at  the  heights 
corresponding  to  the  heights  at  the  midpoint  of  each  pathlength  segment  are 
also  tabulated. 


Several  comments  are  Included  in  the  source  statements  to  annotate  the 
various  steps  in  the  program.  All  symbols  and  terms  used  in  the  computer 
program  statements  are  defined  in  Appendix  D. 
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SYMBOLS  USED  IN  COMPUTER  PROGRAM 


i 


A atmospheric  sound  absorption  coefficient  from  subroutine  ANSAB, 

dB/m 


A1 


All 

A2 

A21 

ACOS 

AE 


atmospheric  sound  absorption  coefficient  along  one  segment  of  a 
propagation  path,  for  test  or  measured  meteorological  conditions, 
and  a frequency  equal  to  f ^ , dB/m 

same  as  Al,  but  for  a frequency  equal  to  dB/m 

same  as  Al,  but  for  reference  meteorological  conditions,  dB/m 

same  as  A2,  but  for  a frequency  of  dB/m 

the  FORTRAN  library  function  arcoslne 

the  exponent  for  the  absorption  term  A(f .)  defined  by  equation 
(B28} , dimensionless  ^ 


AH 

AHM 

AHR 

AJ 


airplane  height  above  ground  level  at  the  time  when  the  sound 
signal  was  emitted  from  the  airplane  noise  source,  m 

height  values  which  are  midway  between  the  heights  for  the 
measurements  of  meteorological  data,  m 


average  relative  humidity  at  the  heights  above  ground  level 
corresponding  to  the  heights  in  the  array  HMI,  percent 


the  absorption  term  A(f  ) in  equation  (B28)  and  equal  to  10' 
dimensionless  ^ 


.A£ 


ALPHA  the  atmospheric  sound  absorption  coefficient  from  equation  (A2) 

and  subroutine  ANSAB,  nepers /m 

AM  airplane  Mach  number  at  the  time  when  the  sound  was  emitted  from 

the  airplane  noise  source,  dimensionless 

AMH  height,  or  vertical  distance,  between  the  airplane  and  the 

microphone,  m 

ANSAB  name  of  a subroutine  for  calculating  pure-tone  sound  absorption 

coefficients  given  the  temperature  (in  kelvlns) , relative  humidity 
(in  percent),  and  pressure  (in  standard  atmospheres)  of  the  air 
and  the  frequency  (in  hertz)  of  the  sound  wave 

AS  airspeed  of  the  airplane  at  the  time  when  the  sound  was  emitted 

from  the  airplane  noise  source,  m/s 
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ATK 

BLJl 

or 

BL3J 

BLTR 

BP 

Cl,  C2, 
C3,  C4 

CF 

D 

DL 

.EQ. 

EXP 

F 

FI 

FJ 

FJl 

FJI 

FRN2 


average  air  temperature  at  the  height  above  ground  level 
corresponding  to  heights  in  the  array  HMI,  kelvins 

band  loss  adjustment  factor  for  the  j-th  frequency  segment  or 
the  sum  over  the  SE  frequency  segments  and  calculated  using 
equations  (B66) , (B67) , or  (B65) , dimensionless 

band  loss  adjustment  for  differences  in  atmospheric  absorption 
between  test  and  reference  conditions,  dB 

measured  station  barometric  pressure  at  the  elevation  of  the 
runway  where  the  noise  tests  were  conducted,  in.  of  Hg  or  mb 
or  kPa 


abbreviations  for  nondimens ional  expressions  used  in  calculating 
the  angle  of  sound  propagation  between  the  aircraft  flight  path 
and  a sound  ray  from  the  aircraft  noise  source  and  the  microphone 

nominal  frequencies  used  for  designating  the  center  frequencies 
of  standard  1/ 3-octave-band  filters,  Hz 


lengths  of  the  segments  of  the  sound  propagation  path  defined  i 

by  the  heights  of  the  measurements  of  meteorological  data  (called  | 

A?  in  equations  (B28)  and  (B70),  m 

differences  in  noise  levels  across  a 1/ 3-octave-band  from  equations 
(B60),  (B61),  or  (B62) , dB 

FORTRAN  library  function  for  'is  equal  to' 

FORTRAN  library  function  for  the  base  e of  natural  logarithms 
and  used  to  calculate  FRN2  in  subroutine  ANSAB 

frequency  of  a sound  wave  in  subroutine  ANSAB,  Hz 

precise  center  frequency  for  a 1/ 3-octave-band  filter  by 
equation  (B48) , Hz 

the  frequency  at  any  one  of  the  SE  steps  in  the  passband  of  a 
filter  between  f and  f and  calculated  according  to  equation 
(B51),  Hz  ^ 


the  next  higher  increment  in  frequency,  f . ^ , after  f , from 
equation  (B47),  Hz  J J 


the  ratio  of  FJ  to  FI  (f^/f^)  from  equation  (B51),  dimensionless 

the  vibrational  relaxation  frequency  of  nitrogen  molecules  from 
equation  (A4)  and  subroutine  ANSAB,  Hz 
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FR02  the  vibrational  relaxation  frequency  of  oxygen  molecules  from 

equation  (A3)  and  subroutine  ANSAB,  Hz 

•GT.  FORTRAN  library  function  for  'is  greater  than' 

H molar  concentration  of  water  vapor  in  a sample  of  moist  air  from 

equation  (A5)  and  subroutine  ANSAB,  percent 


HM 

HMI 

HR 

HRR 

HRT 

lA 

lAl 

IMAX 

ISBN 

K 

K1 

K2 

L 

LOG 

LOGIO 


height  above  ground  level  at  which  upper-air  meteorological  data 
were  measured,  m 

an  array  of  heights  above  ground  level  that  are  midway  between 
the  heights  at  which  upper-air  meteorological  data  were  measured 
and  equal  to  the  midpoint  heights  in  AHM,  m 

relative  humidity  of  the  atmosphere  measured  at  heights  HM,  percent 

relative  humidity  of  the  atmosphere  under  reference  meteorological 
conditions  at  the  midpoint  heights  in  HMI  and  equal  to  the  reference 
relative  humidities  in  the  array  RHR,  percent 

relative  humidity  of  the  atmosphere  under  the  test  (measured) 
meteorological  conditions  at  the  midpoint  heights  in  HMI  and  equal 
to  the  test  relative  humidities  in  the  array  AHR,  percent 

an  index  counter  from  the  array  of  heights  HM  that  indicates  the 
height  for  the  set  of  meteorological  data  that  is  just  less  than 
the  airplane  height  defined  by  AH 

IA-1 

the  number  of  heights  (or  sets  of  meteorological  data)  in  the  array 
HM 

International  Standard  Band  Numbers  for  1/ 3-octave-band  filters 

the  absorption  slope  parameter  over  the  propagation  path  from 
equations  (B29)  and  (B70) , dimensionless 

an  absorption  slope  term  from  equation  (B70)  for  one  segment 
of  the  propagation  path  and  for  the  measured  atmospheric 
conditions,  dimensionless 

same  as  K1  but  for  reference  atmospheric  conditions,  dimensionless 

the  exponent  from  equations  (B60) , (B61),  and  (B62)  for  the  slope 
of  the  measured  1/ 3-octave-band  sound  pressure  levels,  dimensionless 

FORTRAN  library  function  for  natural  logarithms  to  base  e 

FORTRAN  library  function  for  logarithms  to  base  10 


79 


r 


APPENDIX  D 


LPA 

LPSOPO 

.LT. 

.NE. 

PA 

PAR 

PAT 

PD 

PSI 

PSID 

PSOPO 

RBPO 

RF 

RHM 

RHR 

RHRL 

RHRO 

RPA 


the  logarithm  (base  10)  of  the  atmospheric  pressure  (measured 
in  standard  atmospheres  relative  to  1.01325  x lo'  pascals) 
under  reference  meteorological  conditions  at  the  midpoint 
heights  in  HMI,  dimensionless 


as  used  in  subroutine  ANSAB,  the  logarithm  (base  10)  of  the 
ratio  of  the  saturation  vapor  pressure  to  the  reference  atmospheric 
pressure,  dimensionless 

FORTRAN  library  function  for  'is  less  than' 

FORTRAN  library  function  for  'is  not  equal  to' 


measured  air  pressure  in  standard  atmospheres  or  the  ratio  of  BP 
to  an  appropriate  reference  pressure  RBPO,  standard  atmospheres 

air  pressure  in  standard  atmospheres  under  the  reference  meteorological 
conditions  at  the  midpoint  heights  in  HMI  and  equal  to  the  reference 
air  pressures  defined  by  the  array  RPA,  standard  atmospheres 

air  pressure  in  standard  atmospheres  under  the  test  meteorological 
conditions  and  equal  to  PA,  standard  atmospheres 

the  total  length  of  the  sound  propagation  path  as  determined  from 
PD  = AMH/SIN(PSI) , m 


sound  propagation  angle  between  the  airplane's  flight  path  and 
a ray  from  the  source  of  sound  to  the  microphone,  radians 

same  as  PSI  but  in  degrees 


antilog  of  LPSOPO  in  subroutine  ANSAB,  dimensionless 

reference  air  pressure  used  to  determine  PA  from  the  measured  value 
of  pressure  BP  (29.921  in.  Hg  or  1013.25  mb  or  101.325  kPa) 

frequency  ratio  for  1/3-octave  bands  and  equal  to  10^^^^,  dimensionless 

a reference  height  above  ground  level  for  the  reference  values 
of  temperature,  relative  humidity,  and  pressure  and  equal  to  10.0  m 

the  array  of  reference  relative  humidities  at  the  midpoint  heights 
in  array  HMI  from  equation  (2) 

the  lapse  rate  of  the  reference  relative  humidity  with  height 
and  equal  to  -0.0065  percentage  points/m 

the  reference  relative  humidity  at  the  reference  height  RHM  and 
equal  to  70.0  percent 


the  array  of  reference  air  pressures  at  the  midpoint  heights  in  the 
array  HMI  from  equations  (3)  and  (4) 
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RPAL  the  lapse  rate  of  the  logarithm  of  the  reference  air  pressure 

with  helgnt  and  equal  to  -5.393  x 10”^  m~l 

RTK  the  array  of  reference  air  temperatures  at  the  midpoint  heights 

In  array  HMI  from  equation  (1) 

RTKL  the  lapse  rate  of  the  reference  air  temperature  with  height  and 

equal  to  -0.0065  kelvlns/m 

RTKO  the  reference  air  temperature  at  the  reference  height  RHM  and  equal 

to  298.15  kelvlns 

SE  the  number  of  frequency  segments  In  the  passband  of  a filter  and  che 

number  of  steps  In  the  numerical  Integration  of  equation  (B23)  using 
equations  (B68)  and  (B69) , with  the  value  of  20  for  this  program 

SPL  the  measured  1/ 3-octave-band  sound  pressure  levels  In  decibels  re 

20  pPa,  dB 

SPLA  the  1/ 3-octave-band  sound  pressure  levels  after  adjustment  for 

differences  In  atmospheric  absorption  between  test  and  reference 
meteorological  conditions,  dB 

SQRT  FORTRAN  library  function  for  square  root 

TC  air  temperature  measured  at  heights  HM,  degrees  celslus 

TD  the  total  length  of  the  sound  propagation  path  as  determined  by 

summing  the  lengths  of  the  segments  D from  1 to  lA,  m 

TK  air  temperature  used  by  subroutine  ANSAB  in  the  calculation  of 

LPSOPO,  FRN2,  and  ALPHA,  kelvlns 

TKR  air  temperature  under  reference  meteorological  conditions  at  the 

midpoint  heights  In  HMI  and  equal  to  the  air  temperatures  In  the 
array  RTK,  kelvlns 

TKT  air  temperature  under  the  test  (measured)  meteorological  conditions 

at  the  midpoint  heights  in  HMI  and  equal  to  the  air  temperatures 
In  the  array  ATK,  kelvlns 

TO  a reference  temperature  used  by  subroutine  ANSAB  in  the  calculation 

of  FRN2  and  ALPHA  and  equal  to  293.15  kelvlns 

TOl  a reference  temperature  used  by  subroutine  ANSAB  In  the  calculation 

of  LPSOPO  and  equal  to  273.16  kelvlns 

TN  time  of  the  noise  measurements  represented  by  the  array  SPL  and 

relative  to  the  start  time  TS,  seconds 

TOH  time  of  day  on  a 24-hour  clock  (86,400  s)  when  the  aircraft  was 

over  or  closest  to  the  microphone,  seconds 
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time  of  the  midpoint  of  the  noise  data  sample  on  a 24-hour 
(86,400-s)  clock,  relative  to  the  overhead  time  TOH,  seconds 

time  of  day  on  a 24-hour  (86,400-8)  clock  at  the  start  of 
digitizing  the  first  sample  of  data  from  the  original  analog 
magnetic  tape  recording,  seconds 

height  of  the  microphone  above  ground  level,  m 


